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Prologue 



Supersymmetry and the search for new particles 

Elementary particle physics has made enormous progress in the last decades. The 
electroweak and strong interactions of the fundamental building blocks of matter, 
the quarks and leptons, are now described by the so called Standard Model (SM) 
of particle physics. 

The development of the SM was possible due to intensive efforts and successful 
achievements in experiment and theory. On the one side, theoreticians have pro- 
vided the physical models and the mathematical techniques necessary to define 
and calculate observables. On the other side, experiments with particle accelera- 
tors and detectors have not only allowed to find new fundamental particles, but 
also high precision measurements have made it possible to test the models. 

However, there are general arguments which point towards the existence of a the- 
ory beyond the SM. One of the most attractive candidates for such a more funda- 
mental theory is Supersymmetry (SUSY). SUSY transformations change the spin 
of a particle field, and thus bosonic and fermionic degrees of freedom get related 
to each other. Therefore, new particles are predicted in SUSY models like in the 
minimal supersymmetric extension of the SM. If SUSY is realized in nature, the 
supersymmetric partners of the SM particles have to be discovered. 

The next future colliders, like the Large Hadron Collider (LHC) at CERN or a 
planned International Linear Collider (ILC), are designed to find these particles. 
Their properties will be measured with high precision in production and decay 
processes. The underlying physical model can then be determined by a compari- 
son of experimental and phenomenological studies. 



Chapter 1 



Introduction 



1.1 Motivation: Symmetries and models 

Symmetries in physics have always played an important role in understanding 
the structure of the underlying theories. For instance, the existence of conservation 
laws can be explained by specific symmetry transformations under which a theory 
is invariant. Energy momentum and angular momentum are conserved in field 
theories with continuous spacetime symmetries. 

High energy models of elementary particle interactions have to be invariant under 
the transformations of the Poincare group, which are relativistic generalizations of 
the spacetime symmetries. The interactions of the particles, the strong and elec- 
troweak forces, can be understood as a consequence of the so called gauge sym- 
metries. 

But not only such continuous symmetries are crucial. The discrete symmetries 

• Charge conjugation C: interchange of particles with antiparticles 

• Parity P: transformation of the space coordinates x — > — x 

• Time reversal T: transformation of the parameter time t — > —t 

are also essential for the formulation of relativistic quantum field theories. In par- 
ticular, any Poincare invariant local field theory has to be symmetric under the 
combined transformation of C, P and T, which is called CPT invariance. 



1.2 CP violating phases and electric dipole moments 
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The discovery in the fifties that weak interactions violate C and P maximally, was 
noted with substantial belief that a particle theory apparently still conserves the 
combined symmetry CP. However, a few years later, in 1964, CP violation was 
confirmed in the K-meson system. 

This lead to a powerful prediction in 1972. The implementation of CP violation 
made it necessary to include a phase in the quark mixing matrix, the so called 
Kobayashi-Maskawa Matrix [1]. Such a phase is only CP violating if the matrix is 
at least of dimension three, and thus a third generation of quarks and leptons was 
required. This prediction of a third family was made long before the final member 
of the second family, the charm quark, was found. 

Current experiments with B-mesons verify the existence of one CP phase in the 
quark mixing matrix of the SM. However, one phase alone cannot explain the ob- 
served baryon asymmetry of the universe, as shown in [2]. The fact that further 
sources of CP violation are needed leads to the crucial prediction of CP violating 
phases in theories beyond the SM. 

Apart from that, a further important symmetry of physics was born in 1974: Super- 
symmetry (SUSY), which relates fermionic and bosonic degrees of freedom. SUSY 
models, like the the Minimal Supersymmetric Standard Model (MSSM), are one 
of the most attractive theories beyond the SM. They give a natural solution to the 
hierarchy problem and provide neutralinos as dark matter candidates. Further- 
more, SUSY allows for grand unifications and for theories, which might include 
also gravity. 



1.2 CP violating phases and electric dipole moments 

The MSSM might have several complex parameters, which cause CP violating ef- 
fects. In the neutralino and chargino sector these are the Higgsino mass param- 
eter \x = \ii\e lLpti and the U(l) gaugino mass parameter Mi = \Mi\e llfM ^ [3]. The 
SU (2) gaugino mass parameter M 2 can be made real by redefining the fields. In 
the sfermion sector of the MSSM, also the trilinear scalar coupling parameter A f 
of the sfermion / can be complex, A s = \Af\e lpA f . 

The CP violating phases are constrained by electric dipole moments (EDMs) [4] of 
electron e, neutron n, 199 Hg and 205 T1 atoms. Their upper bounds are, respectively: 

|4| < 4.3 x 10~ 27 ecm [5], (1.1) 
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\d n \ < 6.3 x 1(T 26 ecm [6], 
\d Hg \ < 2.1 x 10" 28 ecm [7], 
\d Tl \ < 1.3 x 10- 24 ecm [8]. 



(1.2) 
(1.3) 
(1.4) 



The CP phase in the quark sector of the SM give contributions to EDMs, which 
generally arise at two loop level, and respect the bounds of the EDMs. In SUSY 
models, however, neutralino and chargino contributions to the electron EDM can 
occur at one loop level, see Fig. 1.1. For the neutron EDM in addition also gluino 
exchange contributions are present due to a phase of the gluino mass parameter. 
The phases of the SUSY parameters are thus constrained by the experimental up- 
per limits of the EDMs. In the literature, three solutions are being proposed [9]: 



• The SUSY phases are severely suppressed [10, 11]. 

• SUSY particles of the first two generations are rather heavy, with masses of 
the order of a TeV [12]. 

• There are strong cancellations between the different SUSY contributions to 
the EDMs, allowing a SUSY particle spectrum of the order of a few 100 GeV 
[13-15]. 



Due to such cancellations, for example, in the constrained MSSM [14], the phase 
(fiMi is not restricted but the phase of /i is still constrained with j^j < 0.l7r [14]. 
If lepton flavor violating terms are included [15], also the restriction on y?^ may 
disappear. 

The restrictions on the SUSY phases are thus very model dependent. Independent 
methods for their measurements are desirable, in order to clarify the situation. 
In order to determine the phases unambiguously, measurements of CP sensitive 
observables are necessary. Such observables are non-zero only if CP is violated, 
i.e. they are proportional to the sine of the phases. 




Figure 1.1: SUSY contributions to the electron EDM. 



1.3 Methods for analyzing CP violating phases 
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1.3 Methods for analyzing CP violating phases 
1.3.1 Neutralino and chargino polarizations 

For measuring SUSY phases, the study of neutralino and chargino production at 
an e + e _ linear collider with longitudinally polarized beams [16] will play an im- 
portant role. By measurements of the chargino masses and production cross sec- 
tions, a method has been developed [17-19] to determine cos(y? M ), in addition to the 
other parameters M 2 , and tan/3 of the chargino sector. For neutralino produc- 
tion analogous methods have been proposed in [19-22] to determine also cos(</?mJ 
and Mi, besides cos(y? M ), M 2r \fi\ and tan (3. 

However, in order to determine also the sign of y? M and (fM 1 , the transverse neu- 
tralino and chargino polarizations perpendicular to the production plane have to 
be taken into account [17,18,23]. They are only present if there are CP violat- 
ing phases in the neutralino /chargino sector, and if a pair of different neutrali- 
nos/ charginos is produced. At tree level, their polarizations lead to triple-product 
asymmetries of the decay products [24,25]. Energy distributions and polar angle 
distributions of the neutralinos and charginos or their decay products are not CP 
sensitive at tree level, since they do not depend on the transverse neutralino or 
chargino polarizations, see e.g. [26] for neutralino production. 

In order to include the particle polarizations in our calculations, we use the spin 
density matrix formalism of [27]. For an introduction into this formalism and for 
our conventions and definitions used, see Appendices C and D. 



1.3.2 T-odd and CP-odd triple-product asymmetries 

The SUSY phases give rise to T-odd and CP-odd observables which involve triple 
products of momenta [28]. They allow us to define various T and CP asymmetries 
which are sensitive to the different SUSY phases. On the one hand, these observ- 
ables can be large because they are present at tree level. On the other hand, they 
also allow a determination of the sign of the phases, which is impossible if only 
CP-even observables were studied. 

We consider neutralino or chargino production 

e + + e~ -> Xi + Xj (1-5) 
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followed by the two-body decay of one neutralino or chargino into a SM particle 
A (e.g. lepton or W, Z boson) and a SUSY particle X (e.g. slepton or xT) : 

Xi^A + X. (1.6) 



The momenta of electron, chargino (or neutralino), and particle A define the triple 
product 

T = (p e - xp x J-p A , (1.7) 

which is T-odd, i.e. changes sign under time reversal. The T-odd asymmetry of 
the cross section a of production (1.5) and decay (1.6) is then defined as 

T _ a(T > 0) - ajT < 0) 

~ a(T>0) + a(T < 0) ' 1 j 



The asymmetry can be expressed by the angular distribution of particle A 
J? -^-A cos 9 - Jo 1 -^-A cos 6 N+-N- 

A 1 _ •> i a cos J U acosfl _ + /i n\ 

~ ■ jo jk>dcos9 + r- 1 ^rfcos^ ~ iV+ + AL ' ' 

J i a cost? J U dcosS T 



where 

cosfl:= P«- X P* . -Ed-, (1.10) 

|Pe- X PxJ |PA| 

and thus „4 T is the difference of the number of events with particle A above (JV+) 
and below (AL) the production plane, defined by p e x p Xi , normalized by the 
total number of events iV = N + + AL . 

The T-odd asymmetry is not only sensitive to CP phases, but also to absorptive 
contributions, which could enter via s-channel resonances or final state interac- 
tions at loop level. Although the absorptive contributions are a higher order effect, 
and thus expected to be small, they do not signal CP violation. However, they can 
be eliminated in the CP-odd asymmetry 

-4 CP = ^(-4 T --4 T ), (1.11) 



1.3 Methods for analyzing CP violating phases 
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where A T denotes the asymmetry for the CP-conjugated process. 

Note that the triple product T (1.7) requires the identification of the neutralino (or 
chargino) momentum p Xi , which could be reconstructed by measuring the decay 
of the other Xj- Therefore, the masses of the neutralinos/charginos as well as the 
masses of their decay products have to be known. 

To avoid the reconstruction of p Xi , we can also define triple products in which p x . 
is replaced by a momentum of the decay products of particles A, if A is a W or Z 
boson, or X, if X is a slepton. In the first case, A = W or Z and X = x\, the decay 
of the boson into two quarks 

A^q + q', (1.12) 
defines the triple product 

T = (Pe- xp g )-P q >- (1.13) 
In the second case, A = £ and X — t, the decay 

X^t + xl (1-14) 
defines the triple product 

T = ( Pe - x p A ) ■ p e . (1.15) 

These triple products define then corresponding T or CP asymmetries, which do 
not require the momentum reconstruction of Xi- However, for these triple products 
the leptons have to be distinguished by their energy distributions [29], and the 
quarks have to be distinguished by flavor tagging [30-32]. 

Triple-product asymmetries can also be defined and analyzed for three-body de- 
cays of neutralinos [23,24,33,34] and charginos [25,33,35]. 



1.3.3 Statistical error and significance 



The T-odd and CP-odd asymmetries, as defined in (1.8) and (1.11), could be mea- 
sured in neutralino and chargino production at future linear collider experiments, 
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and would allow us to determine the values of the SUSY phases. In order to de- 
cide whether an asymmetry, and thus a CP phase can be measured, we have to 
calculate its statistical error. Also, we have to consider the statistical significance 
of the asymmetry. 

The relative statistical error of the asymmetry is given by 

SA:=$^ = (1.16) 

1-41 \A\y/N' 

with the number of events N = £ ■ a, where C is the integrated luminosity of the 
linear collider. Formula (1.16) follows from (1.9), with the estimate AN± = y/N± « 

y/N/2. 

The statistical significance of the asymmetry (1.8) is then defined as 



S=\A T \VC^. (1.17) 

For S = 1, the asymmetry can be measured at the 68% confidence level (CL), for 
S = 1.96 at the 95% CL, etc. The significance for the CP-odd asymmetry (1.11) is 
given by 



S=\A CP \y/2C7a, (1.18) 
since A.4 CP = AA T /y/2, which follows from (1.11). 

Also background and detector simulations have to be taken into account to predict 
the expected accuracies for the asymmetries, see e.g. [36]. However, this would im- 
ply detailed Monte Carlo studies, which is beyond the scope of the present work. 



1.4 Organization of the work 

The goal of the thesis is to analyze CP violating effects of MSSM phases in produc- 
tion and/ or two-body decay processes of neutralinos, charginos and sfermions. 
We will therefore define and calculate T-odd and CP-odd asymmetries for the dif- 
ferent supersymmetric processes. 



1.4 Organization of the work 
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We study neutralino and chargino production in electron-positron collisions at a 
future linear collider (LC) with a center of mass energy of 500 GeV to 800 GeV, high 
luminosity and longitudinally polarized beams. A LC of this kind is an ideal tool 
for measuring the properties of SUSY particles with high precision. 

Finally, we address the question, whether the phases can be constrained at the LC. 
We thus calculate the statistical significances for measuring the asymmetries. Our 
analyses will also have particular emphasis on the beam-polarization dependence 
of the asymmetries and cross sections. 

In most of the numerical examples we choose (fM 1 = ±7r/2, (p^ = 0, which is 
allowed by the constraints from the electron and neutron EDMs. In order to show 
the full phase dependences of the asymmetries in some examples we study their 
(p^ behavior in the whole y? M range, relaxing in this case the restrictions from the 
EDMs. This is justified e.g. in theories with lepton flavor violation [15], where the 
constraints on y?^ disappear. 

• Chapter 2 contains neutralino production e + e _ — > Xi X° and decay: 

- In Section 2.1 we discuss neutralino decay into sleptons: xl ~^ ? ? f° r 

£ — e, /i, t. 

- In Section 2.2 we discuss neutralino decay into a stau-tau pair: xl —> r 
including the r polarization. 

- In Section 2.3 we discuss neutralino decay into a Z boson: x° — ► X n 

• Chapter 3 deals with chargino production e + e~ — > xt xj an d decay: 

- In Section 3.1 we discuss chargino decay into sneutrinos: xt —> ^ &e, 

£ — e, /i, t. 

- In Section 3.2 we discuss chargino decay into a W boson: xt ~^ Xn W + . 

• In chapter 4 we analyze CP violation in the two-body decay chain of a sfermion: 

f^f%,7<S^7&Z,Z^ff. 

• Chapter 5 contains a summary and conclusions. 

• In the Appendices we give a short account on the MSSM, with emphasis 
on its complex parameters. We discuss details of particle kinematics and 
phase space and give the analytical formulae for the production and decay 
amplitudes squared. Finally, we give useful spin-formulae for fermions and 
bosons and a formulary of our definitions and conventions, used for our 
numerical calculations. 



Chapter 2 



CP violation in production and decay 

of neutralinos 



Overview 



We study neutralino production with longitudinally polarized beams e + e — > 
X°iX°j with the subsequent leptonic decay of one neutralino xl ~ > W'r £ ~ > Xi^ 
for i = e, r, [29] or the decay into the Z boson x° t -»• X n Z; Z -> [37]. 
These decay modes allow the definition of CP observables which are sensitive to 
the phases y?Mi and y? M . 

For the leptonic decay of the neutralino into the tau Xi —> ^ ± rT > we propose the 
transverse r T polarization as a CP sensitive observable [38]. This asymmetry is 
also sensitive to the phase ^a t - 

We present numerical results for the asymmetries, cross sections and branching 
ratios for a linear electron-positron collider in the 500 GeV - 800 GeV range. The 
asymmetries can go up to 60% and we estimate the event rates which are nec- 
essary to observe the asymmetries. Polarized electron and positron beams can 
significantly enhance the asymmetries and cross sections. 



2.1 T-odd asymmetries in neutralino production and decay into sleptons 
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2.1 T-odd asymmetries in neutralino production and 
decay into sleptons 




Figure 2.1: Schematic picture of the neutralino production and decay process. 
For neutralino production 

(2.1) 



+ e 



with longitudinally polarized beams and the subsequent leptonic two-body decay 
of one of the neutralinos 



Xi - l + ti, 

we introduce the triple-product 

7} = (p e - x p x o) • p £l , 

and define the T-odd asymmetry 

, T = > °) ~ <>A?i < 0) 

1 ar(Tj > 0) + ajiTj < 0) ' 

where 07 is the cross section for reactions (2.1) and (2.2). 



(2.2) 



(2.3) 



(2.4) 
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With the subsequent leptonic decay of the slepton 

€->Xi+£ a; € = e,Ai,T, (2.5) 

we can construct a further asymmetry which does not require the identification of 
the neutralino momentum. We replace the neutralino momentum p x o in (2.3) by 
the momentum p i2 of the lepton from the slepton decay 

T n = (p e - x p h ) ■ p tl (2.6) 
and define the asymmetry 



T = a n (T n > 0) - a n (J n < 0) 
11 a II (T II >0) + a II (T II <0y 

where an is the cross section for reactions (2.1) - (2.5). 



(2.7) 



These T-odd observables in the production of neutralinos at tree level are due to 
spin effects. Only if there are CP- violating phases ip Ml and y? M in the neutralino 
sector and if two different neutralinos are produced, each of them has a polariza- 
tion perpendicular to the production plane [20,23,39]. This polarization leads to 
asymmetries in the angular distributions of the decay products, as defined in (2.4) 
and (2.7). 



2.1.1 Cross sections 



In order to calculate the production and decay amplitudes, we use the spin density 
matrix formalism of [27,39], see Appendix C. For neutralino production (2.1) and 
decay (2.2), the amplitude squared can be written as 

\Ti\ 2 = |A(x°)| 2 £ Pp(x°) AiA Wx°)a^, (2.8) 

with the neutralino propagator A(x° ), the spin-density matrix of neutralino pro- 
duction pp(Xi), the decay matrix pDi(Xi)> an d the neutralino helicities A,, A-. In- 
serting the expansions of the density matrices pp(Xi)> see (CIO), and PdAxD' 
see (C.28), into (2.8) gives 

\Tj\ 2 = 4|A(x°)| 2 (P£i + S P S Dl ). (2.9) 



2.1 T-odd asymmetries in neutralino production and decay into sleptons 
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Analogously, the amplitude squared for the complete process of neutralino pro- 
duction, followed by the two-body decays (2.2) and (2.5), can be written as 

\Tn\ 2 = |A(;a| 2 |A(£)| 2 £ Mx^ <A Wx°)a^ D 2 (2.10) 
= 4 |A(x°)| 2 \A(i)\ 2 (PD 1 + £ P E Dl ) D 2 , (2.11) 
where D 2 is the factor for the slepton decay, given in (C.35). 

The cross sections and distributions in the laboratory system are then obtained by 
integrating the squared amplitudes 

dcriji = ^{Tjj^dUpSjj! (2.12) 
over the Lorentz invariant phase space elements 



dLips 7 := dUps(s;p x o,p ei ,pi), 
dLips /7 := dUps(s;p x o,p £l ,p £2 ,p x o), 

given in (B.22) and (B.23), respectively. 



(2.13) 
(2.14) 



The contributions of the spin correlation terms T, P T lDl to the total cross section 
vanish. Their contributions to the energy distributions of the leptons i\ and i 2 
from decay (2.2) and (2.5) vanish due to the Majorana properties of the neutralinos 
[40] if CP is conserved. In our case of CP violation, they vanish to leading order 
perturbation theory [40], and thus the contributions can be neglected since they 
are proportional to the widths of the exchanged particles. 



2.1.2 T-odd asymmetries 

Inserting the cross sections (2.12) in the definitions of the asymmetries (2.4) and 
(2.7) we obtain 



/Sign[T 7>7/ ]|T| 


2 dUpsjjj 


/ \Ti,ii< 


2 dLips I:II 



aT _ J -o--L-i,iiJ|- I r~l,ll _ J --CD-I- £— uj r~l,ll (r) 1 c\ 
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CP violation in production and decay of neutralinos 



where we have used the narrow width approximation for the propagators. In the 
numerator only the spin-correlation term EpEf^ remains, since only this term con- 
tains the triple products (2.3) or (2.6). Thus, the contributions to Aj n are directly 
proportional to the neutralino polarization Ep perpendicular to the production 
plane. 

In case the neutralino decays into a scalar tau, we take stau mixing into account 
and the asymmetries are reduced due to their dependence on the y^-f^-r couplings 



\ a ki' 


2 - 




2 


14*1 


2 + 




2 



(2.16) 



which can be seen from the expressions of D x and T? Di , given in Appendix (C.2). 
Since the asymmetries are proportional to the absolute values of a T ki , b T ki , they are 
not sensitive to the phase y?A T of A T . As an observable which is sensitive to <pa t , 
we will consider in Section 2.2 an asymmetry which involves the transverse r po- 
larization. 



2.1.3 Numerical results 

We analyze the dependence of the asymmetries Aj and Aj T , the neutralino pro- 
duction cross sections a P (e + e~ — > XiX^j) an d the branching ratios BR(x° — * ££) 
on the parameters /x = e iLp », Mi = \Mi\ e ilfM i and M 2 for tan/3 = 10. In or- 
der to reduce the number of parameters, we assume |Mi] = 5/3M 2 tan 2 6 W and 
fix the universal scalar mass parameter m = 100 GeV. The renormalization group 
equations for the slepton masses are given in (A.56) and (A.57). Since the pair pro- 
duction of equal neutralinos is not CP sensitive, we discuss the lightest pairs x\ x\> 
Xi X3 an d X2 X% f° r which we choose a center of mass energy of y/s — 500 GeV and 
longitudinal beam polarization P e - = 0.8 and P e + = —0.6. 

■ Production ofxi X2 

In Fig. 2.2a we show the cross section for X1X2 production for ip^ — and y?Mi = 
0.5 7r in the \/j,\-M 2 plane. The cross section reaches values up to 300 fb. For 
\fx\ < 250 GeV the right selectron exchange dominates so that our choice of po- 
larizations (P e -,P e +) = (0.8, —0.6) enhances the cross section by a factor up to 2.5 
compared to the unpolarized case. For > 300 GeV the left selectron exchange 
dominates because of the larger X2 ~ &l coupling. In this region a sign reversal of 



2.1 T-odd asymmetries in neutralino production and decay into sleptons 
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both polarizations, i.e. (P e -,P e +) = (—0.8, 0.6), would enhance the cross section 
up to a factor of 20. 

The branching ratio BR(x2 — * £r£i) for the neutralino two-body decay into right 
selectrons and smuons, summed over both signs of charge, is shown in Fig. 2.2b. It 
reaches values up to 64% and decreases with increasing when the two-body de- 
cays into the lightest neutral Higgs boson and/ or the Z boson are kinematically 
allowed. The decays ~^ xfW T are not allowed. With our choice m = 100 GeV, 
the decays into left selectrons and smuons can be neglected because these chan- 
nels are either not open or the branching ratio is smaller than 1%. As we assume 
that the squarks and the other Higgs bosons are heavy, the decay into the stau is 
competing, and dominates for M 2 < 200 GeV in our scenario, see Fig. 2.4a, which 
is discussed below. The resulting cross section is shown in Fig. 2.2c. 

Fig. 2.2d shows the |//|-M 2 dependence of the asymmetry Aj T for ip Ml = 0.5 n and 
y? M = 0. In the region \fj,\ < 250 GeV, where the right selectron exchange dominates, 
the asymmetry reaches 9.5% for our choice of beam polarization. This enhances the 
asymmetry up to a factor of 2 compared to the case of unpolarized beams. With 
increasing the asymmetry decreases and finally changes sign. This is due to 
the increasing contributions of the left selectron exchange which contributes to the 
asymmetry with opposite sign and dominates for > 300 GeV. In this region the 
asymmetry could be enhanced up to a factor 2 by reversing the sign of both beam 
polarizations. 

The sensitivity of the cross section a and the asymmetry Aj T on the CP phases is 
shown by contour plots in the ip^ l -ipM l plane for = 240 GeV and M 2 = 400 GeV 
(Fig. 2.3). In our scenario the variation of the cross section, Fig. 2.3a, is more than 
100%. In addition to the CP sensitive observables, the cross section may serve to 
constrain the phases. For unpolarized beams, the cross section would be reduced 
by a factor 0.4. The asymmetry Ajj (Fig. 2.3b) varies between -8.9% and 8.9%. It 
is remarkable that these maximal values are not necessarily obtained for maximal 
CP phases. In our scenario the asymmetry is much more sensitive to variations of 
the phase ip Ml around 0. The reason is that Ajj is proportional to a product of a 
CP odd (Sp) and a CP even factor (Ef^), see (2.15). The CP odd (CP even) factor 
has as sine-like (cosine-like) dependence on the phases. Thus the maximum of 
Ajj is shifted towards y?Mi = in Fig. 2.3b. On the other hand, the asymmetry is 
rather insensitive to </? M . For unpolarized beams this asymmetry would be reduced 
roughly by a factor 0.33. 

The statistical significance for measuring each asymmetry is given by S = \A T \y/~N 
(1.17), with iV = Ca is the number of events with C the total integrated luminosity. 
We show the contour lines for S = 3 and 5 for Ajj in Fig. 2.3c with £ = 500 fb -1 . 
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Figure 2.2: Contour lines of 2.2a: a>(e + e — > X1X2)/ 2.2b: BR(x° — ► ^idi), i — e,pi, 
2.2c: a P (e+e- - X BR(x° - ^1) X BR(|r - x%) with BR(4 - = 1, 

2.2d: the asymmetry Aj T , in the |//|-M 2 plane for <p Ml = 0.5tc, y? M = 0, tan/5 = 10, 
m = 100 GeV, A T = -250 GeV, = 500 GeV and {P e -,P e +) = (0.8, -0.6). The 
area A (B) is kinematically forbidden by m x o + m x o > y/s (m^ > rn x o). The gray 
area is excluded by m + < 104 GeV. 
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Figure 2.3: Contour lines of 2.3a: a P (eV -> X1X2) x BR(x° -> ^1) x BR(4 -> 
Xi^) with BR(£ fi — > x°£ 2 ) = 1/ 2.3b: the asymmetry Aj If 2.3c: the significance 
5, 2.3d: the asymmetry Aj, in the ip^-ip^h plane for M 2 = 400 GeV, = 240 
GeV, tan/5 = 10, m = 100 GeV, A T = -250 GeV, yfe = 500 GeV and (P e -,P e+ ) = 
(0.8, —0.6). In the gray shaded area of 2.3c we have S < 3. For </? Ml , ^ = we have 
m/ =221 GeV, m Y o = 178 GeV and m Y o = 243 GeV. 
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In Fig. 2.3d we also show the asymmetry Aj which is a factor 2.9 larger than Aj If 
because in Aj r the CP-violating effect from the production is partly washed out 
by the kinematics of the slepton decay. However, for a measurement of Aj the 
reconstruction of the xt momentum is necessary. The asymmetry Aj shows a 
similar dependence on the phases as Ajj because both are due to the non vanishing 
neutralino polarization perpendicular to the production plane. It is interesting to 
note that the asymmetries can be sizable for small values of y? M , which is suggested 
by the EDM constraints, see Section 1.2. 

Next we comment on the neutralino decay into the scalar tau and discuss the main 
differences from the decay into the selectron and smuon. In some regions of the 
parameter space, the decay of the neutralino into the lightest stau fi may dominate 
over that into the right selectron and smuon, and may even be the only decay 
channel. In Fig. 2.4a we show contour lines of the branching ratio BR(x2 ~> ^i r ) 
in the \n\-M 2 plane for A T = -250 GeV, tp Ml = 0.5 tt and w = 0. For M 2 < 200 GeV 
the branching ratio BR(x 2 — > T\t) is larger than 80%. However, due to the mixing 
in the stau sector the asymmetry Aj T , Fig. 2.4b, is reduced compared to that in the 
selectron and smuon channels, see Fig. 2. 2d. The reason is the suppression factor 
(l a Ll 2 — \Ki\ 2 ) / i\ a ki\ 2 + \Ki\ 2 ) (2.16), which may be small or even be zero. 

■ Production of X1X3 

We show in Fig. 2.5a and b contour lines of the cross section <jp(e + e~ — > X1X3) x 
BR(x° -> l R ix) x BR(£ R -> x%) with BR(£ R -> x%) = h and of the asymmetry 
Ajj, respectively. The cross section with polarized beams reaches more than 100 fb, 
which is up to a factor 2.5 larger than for unpolarized beams. The asymmetry 
Ajj, shown in Fig. 2.5b, reaches -9.5%. For unpolarized beams this value would 
be reduced by a factor 0.75. For our choice of parameters the cross section and 
the asymmetry for X1X3 production and decay show a similar dependence on M 2 
and \/i I as for X1X2 production, however, the kinematically allowed regions are 
different. We also studied the y? M dependence of Ajj. For = 0.57r(0.l7r) and 
y?Mi = 0, the maximal values of Ajj in the M 2 -\n\ plane are \Ajj\ < 3%(1%). 

■ Production of X2X3 

The production of the neutralino pair e + e — > X2X3 could make it easier to re- 
construct the production plane because both neutralinos decay. This allows one 
to determine also asymmetry Aj, which is a factor 2-3 larger than Ajj. We dis- 
cuss the decay of the heavier neutralino x°/ which has a larger kinematically al- 
lowed region in the \fi\-M 2 plane than that of m Fig- 2.6 we show the pro- 
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Fig. 2.4a \n\ [GeV] Fig. 2.4b [GeV] 

Figure 2.4: Contour lines of 2.4a: BR(%2 — > t\t) and 2.4b: the asymmetry Aj Ir in 
the |/i|-M 2 plane for tp Ml = 0.5tt, tp^ = 0, A T = —250 GeV, tan/5 = 10, m = 100 
GeV, y/s = 500 GeV and (P e -,P e +) = (0.8, -0.6). The area A (B) is kinematically 
forbidden by m x o+m x o > y/s (m fl > m x o). The gray area is excluded by m x ± < 104 
GeV. 

duction cross section <rp(e + e~ — > X2X3) which reaches 100 fb. The cross section 
a P (e + e- ^ x° 2 x° 3 ) x BR(x° 3 ^ W x BR(£ R - x%) with BR(£ R x°A) = h is 
shown in Fig. 2.6b. The asymmetry Aj T is shown in Fig. 2.6d. As to the y2 M depen- 
dence of Aj, we found that for y? M = 0.57r(0.l7r) and ip Ml = 0, \Aj\ can reach 25% 
(2%) in the |//|-M 2 plane. 

■ Energy distributions of the leptons 

In order to measure the asymmetries Aj (2.4) and *4J 7 (2.7), the two leptons £1 and 
i 2 from the neutralino (2.2) and slepton decay (2.5) have to be distinguished. We 
therefore calculate the energy distributions of the leptons from the first and second 
decay vertex in the laboratory system, see Appendix B.2.4. One can distinguish be- 
tween the two leptons event by event, if their energy distributions do not overlap. 
If their energy distributions do overlap, only those leptons can be distinguished, 
whose energies are not both in the overlapping region. 
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Fig. 2.5a [GeV] Fig. 2.5b [GeV] 



Figure 2.5: Contour lines of 2.5a: a P (e + e- -> x BR(x° -> 44) x BR(4 -> 

with BR(^ -> x^ 2 ) = 1 and £ = e,fi, 2.5b: the asymmetry in the |/i|-M 2 
plane for <^ Ml = 0.5vr, = 0, tan (3 = 10, m = 100 GeV, A T = -250 GeV, = 500 
GeV and (P e -,P e +) = (0.8, —0.6). The area A (B) is kinematically forbidden by 
m x o + m x o > y/s (mj; R > m x o). The gray area is excluded by m x ± < 104 GeV. 

We show in Figs. 2.7a - c different types of energy distributions for lepton l x 
(dashed line), and lepton £ 2 (solid line), i — e,fj,, for e + e~ — > X1X2 an d the sub- 
sequent decays x\ ~~ > ^1 an d i — ► xi^2- The parameters tan (3 = 10, M 2 = 300 GeV, 
W = 0, <^Mi = 0.57T, and for = 200, 300 and 500 GeV, are chosen such that the 
slepton mass = 180 GeV is constant, the LSP mass m x o = 140, 145, 150 GeV is 
almost constant whereas the neutralino mass m v o = 185, 240, 300 GeV is increas- 

ing. The mass difference between £ R and Xi decreases (Am = 40, 35, 30 GeV), 
whereas the mass difference between X2 an d £r increases (Am = 5, 60, 120 GeV). 
The endpoints of the energy distributions of the decay leptons depend on these 
mass differences. Thus, in Fig. 2.7a, the second lepton is more energetic than the 
first lepton. The energy distributions do not overlap and thus the two leptons can 
be distinguished by measuring their energies. This also holds for Fig. 2.7c, where 
the first lepton is more energetic than the second one. In Fig. 2.7b the two distribu- 
tions overlap because the mass differences between x\> £r an d x\ are similar. One 
has to apply cuts to distinguish the leptons, which reduce the number of events. 
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2.1.4 Summary of Section 2.1 

We have considered two triple-product asymmetries in neutralino production e + e~ 
X°iX°j an d the subsequent leptonic two-body decay chain of one neutralino x° — > it, 
i —> Xi^ f° r ^ = e, /x, r. These asymmetries are present already at tree level and are 
due to spin effects in the production and decay process of two different neutrali- 
nos. The asymmetries are sensitive to CP-violating phases of the gaugino and 
Higgsino mass parameters M x and/ or /j, in the neutralino production process. 

For the process e + e~ — > X1X2 an d neutralino decay into a right slepton x% ^r^> 
we have shown that the asymmetries can be as large as 25%. They can be enhanced 
using polarized beams, and can be sizable even for a small phases, y? M , Lp Ml ~ O.lir, 
which is suggested by the experimental limits on EDMs. 
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Fig. 2.6c [GeV] Fig. 2.6d \fj,\ [GeV] 



Figure 2.6: Contour lines of 2.6a: o>(e + e — > X2X3)/ 2.6b: a>(e + e — > X2X3) x 
BR(x° -> x BR(4 -> x?^ 2 ) for £ = e,/i, and BR(4 -> %%) = h 2.6c: 

the asymmetry Aj, 2.6d: the asymmetry Aj Ir in the |/i|-M 2 plane for ip Ml = O.bn, 
^ = 0, tan /3 = 10, m = 100 GeV, A T = -250 GeV, y/s = 500 GeV and {P e - , P e +) = 
(0.8, —0.6). The area A (B) is kinematically forbidden by m x o + m x o > y/s (m^ > 
m v o). The gray area is excluded by m + < 104 GeV. 
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Figure 2.7: Energy distributions in the laboratory system for £ 1 (dashed line) and 
i 2 (solid line) for e + e~ — > X1X2 an d the subsequent decays x°2 ~^ ^R^i an d ?r — > 
X%, for M 2 = 300 GeV, m~ tR = 180 GeV, tan/3 = 10 and {|//|, m x? , m x o}/GeV = 
{200, 140, 185}, {300, 145, 240}, {500, 150, 300} in a, b, c respectively. 
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2.2 CP asymmetry in neutralino production and decay 
into polarized taus 

For the two-body decays of neutralinos into sleptons, where the lepton polariza- 
tions are summed, we have shown in the last section that the asymmetries have 
only CP-odd contributions from the neutralino production process. For neutralino 
decay into a tau, the tau polarization allows to define an asymmetry which has 
also CP-odd contributions from the neutralino decay process. This is particularly 
interesting since an asymmetry can be defined, which is sensitive to the CP phase 
of the trilinear scalar coupling parameter A T . 

We consider neutralino production 

e + + e-^x-+X°p i,j = l,-.-,4, (2.17) 
and the subsequent two-body decay of one neutralino into a tau 

X°^t± + t*; m = l,2. (2.18) 
The r polarization is given by [41] 
„ Tr(p<r ) 



Tr(p) 



(2.19) 



with p the hermitean spin density matrix of the t and o; t the Pauli matrices. The 
component P 3 of the polarization vector P = (Pi, P 2 , P3) is the longitudinal po- 
larization, Pi is the transverse polarization in the plane and P 2 is the transverse 
polarization perpendicular to the plane defined by the momenta p T and p e . The 
transverse polarization P 2 is proportional to the triple product 

T T = s T -(p T x Pe -), (2.20) 

where s T is the t spin 3-vector. For its definition in the laboratory system, see 
B.21. In order to eliminate absorptive phases, we define the CP-odd asymmetry 



A CP = ^(P 2 -P 2 ), 



(2.21) 
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where P denotes the r+ polarization in the charge conjugated process x° —> T m r + . 
The asymmetry Aqp is sensitive to the phase <^a t in the stau sector, as well as to 
the phases and y? Ml in the neutralino sector. 

Note, that for the two-body decay (2.18), the transverse r polarization P 2 is the 
only observable which is sensitive to ^ At . The asymmetries defined in Section 2.1, 
where the r polarization is summed, are only sensitive to CP violation due to y? M 
and (p Ml in the production process. 

2.2.1 Tau spin-density matrix and cross section 

In the spin density matrix formalism, see Appendix C, the unnormalized spin- 
density matrix of the r can be written as 

p P (r-)^ = |A(x°)| 2 E PP^f lK PDi£j$. (2-22) 

It is composed of the neutralino propagator A(%°), the spin density matrices pp(Xi) 
for neutralino production (2.17) and Pd(x1.) f° r neutralino decay (2.18). The Xi 
helicities are denoted by Aj and \' ir and the t" helicities are denoted by and X' k . 
The neutralino production matrix pp(x°) is defined in (C.10) and the neutralino 
decay matrix Pd{xD in (C.38). Inserting these density matrices into (2.22) gives 

p P {r-)^ = 4|A(x°)| 2 l(PD + ^ a D )5x kK + (PD b + S« E° b )< A , )]. (2.23) 

The last term of the coefficient S^", see (C.44), contains for 6 = 2 the triple prod- 
uct (2.20). This term is proportional to the product of the x°-r fc -r couplings lm.(b T m * a T mi ) 
and is therefore sensitive to the phases <^> At , y? M and Lp Ml ■ 

The amplitude squared is obtained by summing over the r helicities in (2.23) 
|T| 2 = 4|A(x°)| 2 [P (2D) + E° (2EJ,)], (2.24) 
where the CP sensitive term E^ drops out. The cross section is then given by 

da = ^-\T\ 2 dLips(s;p x o,p T ,p f ), (2.25) 
with the phase space element rfLips as defined in (B.22). 
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2.2.2 Transverse tau polarization and CP asymmetry 

From (2.23) we obtain for the transverse r~ polarization (2.19) 

_ /SfrSgrfLips 
2 JPDdLips ' 



(2.26) 



which follows since we have used the narrow width approximations for the prop- 
agators and in the numerator / |A(x°)| 2 PD 2 rfLips = and in the denominator 

/|A(x°)| 2 S^rfLips = 0. 

As can be seen from (2.26), P 2 is proportional to the spin correlation term £^ 2 (C.44), 
which contains the CP-sensitive part lm(b T m * a T mi )e^ vp(T p^p v - s a \f s b T ' a . In order to 
study the dependence of P 2 on the parameters, we expand for f ± 



Im(6[/atJ = /cos^fYrlm^Ay + g 2 sin 2 ^y T v / 2tan^Im(^ 1 ^ 3 ) 
+ <? 2 sin 2 f cos 2 ^[y T 2 Im(7V i3 ^3e^) + g 2 V2t a n6 w Im(f T Ll N ll e- t ^)}, (2.27) 

using the definitions of the couplings in the stau sector, see Appendix A.2. 

For ipfxjifiMx = 0, we find from (2.27) that P 2 oc sin 26 f sin tp? . We note that the 
dependence of (pf on (p Ar is weak if \A T \ <C tan/3, see (A.49). Thus, we expect 
that P 2 increases with increasing \A T \. 

In order to measure P 2 and the CP asymmetry Aqp (2.21), the t from the neu- 
tralino decay x° — > r+r~ and the r + from the subsequent f+ decay f + — > x? r+ 
have to be distinguished. This can be accomplished by their different energy dis- 
tributions, see Appendix B.2.4. 



2.2.3 Numerical results 

We present numerical results for e + e~ — > X1X2 an d the subsequent decay of the 
neutralino into the lightest stau x° ~~ * f° r a linear collider with y/s = 500 GeV 
and longitudinally polarized beams with (P e -,P e +) = (±0.8,^0.6). This choice 
favors right or left selectron exchange in the neutralino production process, re- 
spectively. 

We study the dependence of the asymmetry .Acp and the production cross sections 

a = a P (e + e~ -> X1X2) x BR(x§ -> f^r - ) on the parameters </v \fi\, ip Ml , \Mi\, <p At , 
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Figure 2.8: Contour lines of Acp for 
M 2 = 200 GeV, = 250 GeV, tan (3 = 
5, (pMi = W = and (P e -,P e +) = 
(0.8,-0.6). 
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Figure 2.9: Contour lines of .Acp 
for A T = 1 TeV, M 2 = 300 GeV, 
\fi\ = 250 GeV, tp Ar — tfp — and 
(P e -,P e+ ) = (0.8,-0.6). 



|A r | and tan/?. We assume jM^ = 5/3 M 2 ta.n 2 6 w , and use m = 100 GeV for 
the universal scalar mass parameter in the renormalization group equations of the 
selectron masses, see (A.46) and (A.47). We take into account the restrictions on 
\A T \ due to the tree-level vacuum stability conditions [42]. 

For the calculation of the branching ratio BR(x° —> t~i t ^) w e include the two-body 
decays 

Xl - r m r, £ R , L £, x°Z, ifiW*, x\Hl £ = e,v, m, n = 1, 2, (2.28) 

with rriA = 1 TeV, such that the neutralino decays into the charged Higgs bosons 
X2 — > xtH T , as well as decays into the heavy neutral Higgs bosons x°2 —> Xi ^2,3/ 
are excluded in our scenarios. 

In Fig. 2.8 we show the contour lines of Acp in the <fA T -\ A T \ plane. The asymmetry 
^4cp is proportional to sin 26 f sin </? f , and increases with increasing \A T \ 3> \fi\ tan (3, 
which is expected from (2.27). Furthermore, in the parameter region shown the 
cross section a varies between 20 fb and 30 fb. 
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Figure 2.10: Contour lines of a and Aqp in the \fx\-M 2 plane for (fA T = 0.57T, Lp Ml = 
y? M = 0, A T = 1 TeV, tan/? = 5 and (P e - ) P e +) = (—0.8, 0.6). The blank area outside 
the area of the contour lines is kinematically forbidden since here either ^/s < 
m x o + m x o or m fl + m T > m x o. The gray area is excluded by m x ± < 104 GeV. 



In Fig. 2.9 we show the dependence of Acp on tan (5 and ip Ml - Large values up to 
±20% are obtained for tan (3 w 5. Note that these values are obtained for 
±0.87r rather than for maximal v?Mi ~ ±0.57r. This is due to the interplay of CP- 
even and CP-odd contributions to the spin correlation terms in (2.26). In the region 
shown in Fig. 2.9, the cross section a varies between 10 fb and 30 fb. 

Figs. 2.10a and 2.10b show, for (p Ar = 0.5n and (p Ml = V 9 ^ — 0/ t ne \fi\-M 2 depen- 
dence of the cross section a and the asymmetry .Agp, respectively. The asymmetry 
reaches values up to —15% due to the large value of \A T \ = 1 TeV and the choice of 
the beam polarization (P e -,P e +) = (—0.8, 0.6). This choice also enhances the cross 
section, which reaches values of more than 100 fb. 

For (p Ml = 0.57T and y? M = ^a t = we show in Figs. 2.11a,b the contour lines of a 
and Acp, respectively, in the \fi\-M 2 plane. It is remarkable that in a large region 
the asymmetry is larger than —10% and reaches values up to —40%. Unpolarized 
beams would reduce „4cp only marginally, however the largest values of a would 
be reduced by a factor 3. 

For \fi\ = 300 GeV and M 2 = 400 GeV, we show in Figs. 2.12a,b contour lines of a 
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Figure 2.11: Contour lines of a and Acp in the \/jl\-M 2 plane for <p Ml = 0.5n, (p Ar = 
tpp = 0, A T = 250 GeV, tan/? = 5 and (P e -,P e +) = (-0.8,0.6). The blank area 
outside the area of the contour lines is kinematically forbidden since here either 

y/s < m x o + m x o or m fl + m T > m x o. The gray area is excluded by m x ± < 104 GeV. 



and Acp, respectively in the (p p -(p Ml plane. The asymmetry Acp is very sensitive 
to variations of the phases ^>m 1 and y? M . Even for small phases, e.g. </? M , ~ 0.1, 
we have .Acp ~ 15%. 

The polarization of the r can be analyzed through its decay distributions. The sen- 
sitivities for measuring the polarization of the r for the various decay modes are 
given in [43]. The numbers quoted there are for an ideal detector and for longitudi- 
nal r polarization and it is expected that the sensitivities for transversely polarized 
r leptons are somewhat smaller. Combining informations from all r decay modes 
a sensitivity of S — 0.35 [44] has been obtained. Following [43], the relative sta- 
tistical error of P 2 (and of P 2 analogously) can be calculated as 5P 2 = AP2/I-P2I = 
a s /(S\P2\ y/N), for a s standard deviations, and N = aC events for the integrated 
luminosity C and the cross section a = ap(e + e~ — > X1X2) x BFi(X2 ~ ^i +r )• Then 
for Acp (2.21), it follows AAcp = AP 2 /\/2. We show in Fig. 2.13 the contour lines 
of the sensitivity S = \/2 / (\Acp\VN) which is needed to measure Acp at 95% CL 
(a s = 2) for C = 500 fbT 1 , for tfA T = 0.2tt and •^>m 1 = W = 0. In Fig. 2.14 we show 
the contour lines of the sensitivity S for ip Ml = 0.27T and (p^ = (fA T = 0. It is inter- 
esting to note that in a large region in the \ji\-M 2 plane in Figs. 2.13 and 2.14 we 
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Figure 2.12: Contour lines of <r and ^4cp in the Lp p -Lp Ml plane for M 2 = 400 GeV, 
= 300 GeV, tan/? = 5, ^ T = 0, A T = 250 GeV and (P e -,P e +) = (-0.8, 0.6). 



obtain a sensitivity 5 < 0.35, which means that the asymmetries can be measured 
at 95% CL. 



2.2.4 Summary of Section 2.2 

We have defined and analyzed a CP odd asymmetry Aqp of the transverse r po- 
larization in neutralino production e + e~ — > an d subsequent two-body decay 
Xi ~^ Tj*rr T . The asymmetry is sensitive to CP- viola ting phases of the the trilin- 
ear scalar coupling parameter A T and the gaugino and Higgsino mass parameters 
Mi, ii. The asymmetry occurs already at tree level and is due to spin effects in the 
neutralino production and decay process. In a numerical study for e + e~ — > X1X2 
and neutralino decay x\ ~^ ^i tT w e have shown that the asymmetry can be as 
large as 60%. It can be sizable even for small phases of \i and M lr suggested by the 
experimental limits on EDMs. 
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Figure 2.13: Contour lines of S for 

(p Ar = 0.2tt, tp Ml = tpp = 0, A T = 1 TeV, 
tan/5 = 5 and (P e -,P e+ ) = (-0.8,0.6). 
The blank area outside the area of the 
contour lines is kinematically forbid- 
den since here either y/s < m x o + m x o 
or m fl + m T > m x o. The gray area is 
excluded by m x ± < 104 GeV. 
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Figure 2.14: Contour lines of S for 

ip Ml = 0.27T, ip Ar = ip^ = 0, A T = 
250 GeV, tan/? = 5 and (P e -,P e+ ) = 
(—0.8, 0.6). The blank area outside the 
area of the contour lines is kinemati- 
cally forbidden since here either y/s < 
m x o + m x o or + m T > m x o . The gray 
area is excluded by m + < 104 GeV. 
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2.3 T-odd observables in neutralino production and 
decay into a Z boson 



Figure 2.15: Schematic picture of the neutralino production and decay process. 

A further possibility to study CP violation in the neutralino sector is the two-body 
decay of the neutralino into a Z boson. Due to spin correlations of the neutralino 
and the Z boson, observables can be defined which have not only CP-odd contri- 
butions from the neutralino production, but also from its decay. 

We study CP violation in neutralino production 




Pf 



e + + e 



Xi+Xj, z,j = l,...,4, 



(2.29) 



with the subsequent two-body decay of one neutralino into a Z boson 



Xi^Xl + Z; n<i, 



(2.30) 



followed by the decay of the Z boson 



Z->f + f] f = t,Q, £ = e,fi,T, q = c,b. 



(2.31) 



For a schematic picture of the neutralino production and decay process see Fig. 2.15. 
If CP is violated, the phases <p Ml and ip^ lead to CP sensitive elements of the Z 
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boson density matrix. They involve CP-odd asymmetries Aj in the angular distri- 
bution of the decay fermions 

_ *(T f > 0) - v(T f < 0) 

Af ~ o(T f >o) + o(T f <oy {I - 3I) 

of the triple product 

T f = p e - • ( P/ x P/ ), (2.33) 

and the cross section a of neutralino production (2.29) and subsequent decay chain 
(2.30)-(2.31). Due to the correlations between the x° polarization and the Z boson 
polarization, there are CP-odd contributions to the Z boson density matrix and to 
the asymmetries Af from the production (2.29) and from the decay process (2.30). 

In Section (2.1) we have studied asymmetries for neutralino decay into sleptons 
X°i — > it. We have shown that these asymmetries have in contrast only contri- 
butions from the neutralino production process, since the neutralino decay is a 
two-body decay into scalars. 

Note that if we would replace the triple product 7} by Tf — p e • (p x <> x p z ), and 
would calculate the corresponding asymmetry, where the Z boson polarization is 
summed, all spin correlations and thus this asymmetry would vanish identically 
because of the Majorana properties of the neutralinos. 



2.3.1 Cross section 



For the calculation of the cross section for the combined process of neutralino pro- 
duction (2.29) and the subsequent two-body decays (2.30), (2.31) of x° we use 
the same spin-density matrix formalism as in [27,39]. The (unnormalized) spin- 
density matrix of the Z boson 

PP (Z)^ = |A(x?)| 2 E Pp^ K PnAx^t (2-34) 



is composed of the neutralino propagator A(y°), the spin-density production ma- 
trix pp{x°i), defined in (C.3), and the decay matrix poAxDr defined in (C.46). The 
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amplitude squared for the complete process e + e — > x°x°; Xi —> x\Z) Z —> ff can 
now be written 

|T| 2 = |A(Z)| 2 £ Pp(Z) x ^ p D2 (Z) x , Xk , (2.35) 

with the decay matrix Pd 2 (Z) for the Z decay, defined in (C.49). Inserting the 
density matrices pp(Xi) (C.10) and Pdi(x°) (C.63) into (2.34) leads to 

p P (Z) XkX * = 4 |A(x-)| 2 [PD X 5 XkX * + Y?p c Y? Di (J c ) AfcA * + P cd D 1 (J"*)****] , (2.36) 

summed over a, c, d. Here the Z production matrix pp(Z) is decomposed into con- 
tributions of scalar (first term), vector (second term) and tensor parts (third term). 
Inserting then p P (Z) (2.36) and p D2 (Z) (C.64) into (2.35) leads finally to 

|T| 2 = 4\A{x^)\ 2 \A(Z)\ 2 [3PD l D 2 + 2Z P c i:% i c D 2 + 4P( cd Dl d D 2 -l cc D 1 dd D 2 

(2.37) 

The differential cross section is then given by 

do= ^-|T| 2 rfLips(s;p x o,p x o,p / ,p / ), (2.38) 
where rfLips is the Lorentz invariant phase-space element defined in (B.40). 



2.3.2 Z boson polarization 

The mean polarization of the Z boson is given by its 3 x 3 density matrix < p(Z) > 
with Tr{< p[Z) >} = 1. We obtain < p(Z) > in the laboratory system by inte- 
grating (2.36) over the Lorentz invariant phase-space element dLips(s; p x q,p x o ,p z ) 
(B.39), and normalizing by the trace 

< P^ kK >= r if ( fiXvwT S = ¥ XkK + V. (J c ) XkK + T cd (J cd ) XkK - (2-39) 
J Tr{pp(z) AfcA fc| rfLips 
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The components V c of the vector polarization and T cd of the tensor polarization are 
given by 

J^y^dUps _ [P cd Di dLips 

Vc - 3 / 2 PD, dLips ' cd ~ dc ~ 3 / PD 1 dLips ' 1 } 

where we have used the narrow width approximation for the neutralino propaga- 
tor. The tensor components T 12 and T 23 vanish due to phase-space integration. The 
density matrix in the helicity basis, see Appendix F.2, is given by 



<p{Z)- > = ±-V 3 + T 33 , (2.41) 

< p (Zf° > = -2T 33 , (2.42) 

< p(Z)-° > = ^-(V 1 + iV 2 ) - V2T 13 , (2.43) 

< p(Z)~ + > = T u , (2.44) 

< p(Z) 0+ > = i(K + iV 2 ) + V2T 13 , (2.45) 



where we have used T u + T 22 + T 33 = -\ and T 12 = T 23 = 0. 



2.3.3 T-odd asymmetry 

From (2.37) we obtain for the asymmetry (2.32) 

_ / Signal 2 ^Lips _ /Signp}] 2 S^S^ C D 2 rfLips 
f /|T| 2 dLips J 3 PD ± D 2 dLips ' l ' ; 

and dLips(s;p x q,p x o,pf,pj) is the Lorentz invariant phase-space element defined 
in (B.40), where we have already used the narrow width approximation for the 
propagators. In the numerator only the vector part of |T| 2 remains which contains 
the triple product 7} = p e ■ (p/ x pj). In the denominator the vector and ten- 
sor parts of |T| 2 vanish, since the complete phase-space integration eliminates the 
spin correlations. Due to the correlations Sp c Ti a Dl between the and the Z boson 
polarization, there are CP-odd contributions to the asymmetry Af from both the 
neutralino production process (2.29), and from the neutralino decay process (2.30). 
The contribution from the production is given by the term with a = 2 in (2.46) and 
is proportional to the transverse polarization Sp (C.17) of the neutralino perpen- 
dicular to the production plane. For the production of a pair of equal neutralinos, 
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e + e — »■ XiXif we have £p = 0. The contributions from the decay, given by the 
terms with a = 1, 3 in (2.46), are proportional to 

C E^ C D 2 d -8m x o (ImO^)(ReO'^)(R 2 - L 2 )(t c z • W<fPx#/, (2-47) 

see last term of (C.62), which contains the e-tensor. Thus Af may be enhanced 
(reduced) if the contributions from production and decay have the same (opposite) 
sign. 

Note that the contributions from the decay would vanish for a two-body decay of 
the neutralino into a scalar particle, as discussed in Section (2.1). We have found 
in this case, that only contributions to Af from the production remain, which are 
multiplied by a decay factor oc (\R\ 2 -\L\ 2 ), and thus A f oc (\R\ 2 -\L\ 2 )/(\R\ 2 +\L\ 2 ), 
see (2.16), where R and L are the right and left couplings of the scalar particle to 
the neutralino. 

For the measurement of Af the charges and the flavors of / and / have to be 
distinguished. For / = e, \x this will be possible on an event by event basis. For / = 
t it will be possible after taking into account corrections due to the reconstruction 
of the r momentum. For / = q the distinction of the quark flavors should be 
possible by flavor tagging in the case q = b,c [30,31]. However, in this case the 
quark charges will be distinguished statistically for a given event sample only. 
Note that A q is always larger than Ai, due to the dependence of Af on the Z-f-f 
couplings, which follows from (C.65), (C.66) and (2.46): 

* -V> = % + f 1 T Ife * * 6-3 (4.5) x At, (2.48) 

U f + L f n t~ L t n b(c) + L b(c) 

compare also Section 4.2. 

The significance for measuring the asymmetry is given by Sf = \Af\y/~N, see (1.17). 
Note that Sf is larger for f = b,c than for / = £ — e, /i, r with Sb ^ 7.7 x Si and 
S c ~ 4.9 x S t , which follows from (2.48) and from BR(Z -> 66) ~ 1.5 x BR(Z -> tl), 
BR(Z -> cc) ~ 1.2 x BR(Z -> [45]. 



2.3.4 Numerical results 

We study the dependence of the Z density matrix < p(Z) > (2.39), the asymme- 
try At(i = e,n,r) (2.32), and the cross section a = a P (e + e~ — > xlx'j) x BR(x° — > 
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jftZ) x BR(Z -> it) on the MSSM parameters \i = e*^ and M 1 = \M ± \ e^ M i 
for tan (3 = 10. In order to reduce the number of parameters, we assume the re- 
lation | Mi | = 5/3M 2 tan 2 9 W and use the renormalization group equations for the 
sfermion masses, Appendix A.3.4, with m = 300 GeV. For the branching ratio 
Z — > it, summed over t = e, fx, r, we take BR(Z — * it) = 0.1 [45]. The values for 
Ab, c are given by (2.48). We choose a center of mass energy of y/s — 800 GeV and 
longitudinally polarized beams with beam polarizations (P e -,P e +) = (±0.8, =f0.6). 

For the calculation of the neutralino widths T x o and branching ratios BR(x° — > 
X?z"), see Appendix E.l, we include the following two-body decays, if kinemati- 
cally allowed, 

Xi -»• g i?,Le, T"mr, i/^, Xm W± » xM, £ = e,fi,r, m = l,2, n< <2.49) 

and neglect three-body decays. The Higgs parameter is chosen = 1 TeV and in 
the stau sector, we fix the trilinear scalar coupling parameter A T = 250 GeV. 

■ Production of Xi X 2 

In Fig. 2.16a we show the cross section for X1X2 production in the \/j,\-M 2 plane for 
ip^ = and ip Mi = 0.57T. For |//| > 250 GeV the left selectron exchange dominates 
due to the larger x 2 — &l coupling, such that the polarization (P e - , P e +) = (—0.8, 0.6) 
enhances the cross section to values of more than 110 fb. The branching ratio 
BR(x° — ► Zxi), see Fig. 2.16b, can even be 100% and decreases with increasing 
|//| and M 2 , when the two-body decays into sleptons and/ or into the lightest neu- 
tral Higgs boson are kinematically allowed. The cross section a = a P (e + e~ — > 
X1X2) x BR(x° — >■ Zxi) x BR(Z — > it), see Fig. 2.16c, does however not exceed 
7 fb, due to the small BR(Z -> It) = 0.1. Fig. 2.16d shows the |/x|-M 2 dependence 
of the asymmetry ^ for (p Ml = 0.5tt and = 0. The asymmetry |^| can reach 
a value of 1.6%. The (positive) contributions from the production cancel the (neg- 
ative) contributions from the decay on the contour At — 0. We also studied the 
ifn dependence of At. In the \\i\-M 2 plane for ip Ml = and p^ = 0.5n we found 
\At\ < 0.5%. 

In Fig. 2.17 we show the ip^-ip Ml dependence of At for |//| = 400 GeV and M 2 = 
250 GeV. The asymmetry At is more sensitive to p Ml than to y? M . It is remarkable 
that the maximal phases of p Ml ,p fl = ±n/2 do not lead to the highest values of 
At ~ ±1.4%, which are reached for ((p Ml ,P^) ~ (±0.37r, 0). The reason for this is 
that the spin-correlation terms £ P c SJ )i C D 2 in the numerator of Af (2.46), are prod- 
ucts of CP-odd and CP-even factors. The CP-odd (CP-even) factors have a sine-like 
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Figure 2.16: Contour plots for 2.16a: a P (e + e- -> X1X2)/ 2.16b: BR(x° -> Z*?), 
2.16c: a = a P (e + e - -> x B R(x° -> ^Xi) x BR(Z -> with BR(Z -> 

££) = 0.1, 2.16d: the asymmetry At, in the |/x|-M 2 plane for cp Ml = 0.5ir, cp^ = 0, 
tan/? = 10, m = 300 GeV, = 800 GeV and (P e -,P e +) = (-0.8, 0.6). The area A 
(B) is kinematically forbidden by m^o + m^o > y/s (mz + m^o > m$o). In area C of 
plot 2.16b: BR(x° -> Zx°) = 100%. The gray area is excluded by m~± < 104 GeV. 
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Figure 2.18: Vector (Vj) and tensor 
(T u ) components of the Z density 
matrix for e+e" -> X1X2, xl ~* z x\, 
for M 2 = 250 GeV, = 400 GeV, 

W = 0, tan (5 = 10, m = 300 
GeV, = 800 GeV and (P e -,P e +) = 
(-0.8,0.6). 



(cosine-like) phase dependence. Therefore, the maximum of the CP asymmetry Af 
is shifted from tp Ml , <fn = ±7r/2 to a smaller or larger value. 

In the f^-fMi region shown in Fig. 2.17 also the cross section er = a P (e + e~ — > 
Xlxl) x BR(x° - Zx?) x BR(zT - with BR(x° - Zxi) = 1 and BR(zT - 
££) = 0.1, is rather insensitive to and varies between 7 fb (</?Mi = 0) and 14 fb 
(y?Mi = ±7r). The statistical significance for measuring the asymmetry for the lep- 
tonic decay of the Z is given by Se = \Ae\VC ■ a, see Section 2.3.3. For C = 500 fb -1 , 
we have Se < 1 in the scenario of Fig. 2.17 and thus At cannot be measured at the 
68% confidence level (Se = 1)- For hadronic decays into b(c) quarks, however, the 
significance is larger S b ^ = 7.7(4.9) as discussed Section 2.3.3. For C = 500 fb -1 
and ((pMi, V/x) — (±0.37r, 0) in Fig. 2.17 we find S b ^ = 8(5) and thus Ab( c ) can be 
measured. 

In Fig. 2.18 we show the dependence of the vector (Vj) (Vj) and tensor (To) 
components of the Z boson polarization. The components Tn, T 22 and V\ have a 
CP-even dependence on (p Ml . The component V 2 is CP-odd and is not only zero 
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for ip Ml = and <pm 1 = n, but also for ip Ml ~ (1 ± 0.2)7r, due to the destructive 
interference of the contributions from CP violation in production and decay. The 
interference of the contributions from the CP-even effects in production and decay 
cause the two maxima of V\. As discussed in Appendix C.4, the tensor components 
Tn and T 22 are almost equal. Compared to V± and V 2 , they have the same order of 
magnitude but their dependence on tp Ml is rather weak. The components Ti 3 , V 3 < 
are small, and thus the density matrix < p{Z) > is almost symmetric. In 
the CP conserving case, e.g. for <p Ml = W = 0, M 2 = 250 GeV, \p\ = 400 GeV, 
tan/5 = 10, m = 300 GeV, = 800 GeV and (P e -,P e +) = (-0.8, 0.6) it reads 

0.329 0.049 0.0003 

< p (Z) >= | 0.049 0.343 0.049 ] . (2.50) 
0.0003 0.049 0.329 

In the CP violating case, e.g. for <p Ml = 0.57T and the other parameters as above, 

< p(Z) > has imaginary parts due to a non-vanishing V 2 



< P{Z) >- 



( 0.324 0.107 + 0.037i 0.0003 \ 

0.107 - 0.037i 0.352 0.107 + 0.037i 
V 0.0003 0.107 - 0.037i 0.324 / 



(2.51) 



Imaginary elements of < p[Z) > are thus an indication of CP violation. Note that 
also the CP even diagonal elements are changed for (p Ml ^ (and also for y? M ^ 0). 
This fact has been exploited in [46] as a possibility to determine the CP violating 
phases. The (p Ml , W dependence of the Z-density matrix elements has also been 
studied in [47], for e + e~ — > X1X3 followed by xl %x\- 



■ Production of \\x\ 



In Fig. 2.19a we show the cross section a = a P (e + e~ — > X2X2) x BR(x 2 ZxT) x 
BR(Z — > it) in the \p\-M 2 plane for ip^ = and v?Mi = 0.57T. The production 
cross section ap(e + e~ — > X2X2)/ which is not shown, is enhanced by the choice 
(P e -,P e +) = (—0.8,0.6) and reaches values up to 130 fb. The branching ratio 
BR(x[] -»• Zxl) can be 100%, see Fig. 2.16b, however, due to BR(Z -> if) = 0.1, 
a is not larger than 13 fb, see Fig. 2.19a. 

If two equal neutralinos are produced, the CP sensitive transverse polarization of 
the neutralinos perpendicular to the production plane vanishes, T? P = in (2.46). 
However, the asymmetry Af obtains CP sensitive contributions from the neu- 
tralino decay process, terms with a = 1,3 in (2.47). In Fig. 2.19b we show for 
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Figure 2.19: Contour lines of a = ap(e + e~ -> xlxl) x B R(x° -» x BR(Z -> 
(2.19a), and the asymmetry ^ (2.19b) in the \/i\-M 2 plane for ip Ml = 0.5tc, tp^ = 0, 
tan p = 10, m = 300 GeV, y/s = 800 GeV and {P e -, P e +) = (-0.8, 0.6). The area A 
(B) is kinematically forbidden by m^o + m^o > y/s {m z + m%o > m^o). 
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Figure 2.20: Vector (Vi) and tensor (Tn) components of the Z density matrix for 

e + e- -> X2X2, X°2 -> Zxl for M 2 = 250 GeV, |//| = 400 GeV, ^ = 0, tan/? = 10, 
m = 300 GeV, = 800 GeV and (P e -,P e+ ) = (-0.8, 0.6). 



V?Mi = 0.57T and tp^ = the and M 2 dependence of the asymmetry At, which 
reaches more than 3%. Along the contour At = in Fig. 2.19b the contribution to 
At which is proportional to Ep, see 2.46, cancels that which is proportional to Ep. 
As the largest values of At > 0.2% and „4 ? > 1% lie in a region of the \fi\-M 2 plane 
where o < 0.3 fb, it will be difficult to measure Af in a statistically significant way. 
We also studied the dependence of In the \/jl\-M 2 plane for (p Ml = and 
p^ = 0.57T we found |^| < 0.5%, and thus the influence of p^ is also small. 



In Fig. 2.20 we show the p>M x dependence of the vector (Vi) and tensor (Tn) com- 
ponents of the Z boson polarization. Because there are only CP sensitive contribu- 
tions from the neutralino decay process, V 2 is only zero at = 0, n and V\ has one 
maximum at (f Ml = n, compared to the components for xi x\ production, shown 
in Fig. 2.18. In addition, the vector components V 1 and V 2 in Fig. 2.20 are much 
smaller than the tensor components Tn fx T 22 . The smallness of V 2 accounts for 
the smallness of the asymmetry \At\ < 0.05%. Furthermore, the other components 
are small, i.e. T 13 < 10" 6 and V 3 = 0. 
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Figure 2.21: Contour lines of a = a P (e + e~ -> xlxi) x BR(x° -> ^X?) x BR(Z -> ^) 
(2.21a), and the asymmetry Ae (2.21b) in the \fi\-M 2 plane for (p Ml = 0-5tt, (p^ = 0, 
tan/3 = 10, m = 300 GeV, yfi = 800 GeV and (P e -,P e +) = (0.8, -0.6). The area 
A (B) is kinematically forbidden by m^o + m^o > y^s (m z + m^o > m%o). The gray 
area is excluded by m~± < 104 GeV. 



Production of X\ X.3 



xixi) 



x BR(x° -> ZxV 



In Fig. 2.21a we show the cross section a = ap(e + e~ 
BR(Z — > ££) in the |/i|-M 2 plane for = and <£>Mi = 0.57T. The production 
cross section o>(e + e~ — > X1X3)/ which is not shown, is enhanced by the choice 
(P e -,P e +) — (0-8, —0.6) and reaches up to 50 fb. The branching ratio BR(x° — > 
Zxi), which is not shown, can be 100%, however, due to BR(Z — > ££) = 0.1, the 
cross section shown in Fig. 2.21a does not exceed 5 fb. In Fig. 2.21b we show the 
|//|-M 2 dependence of the asymmetry At. The asymmetry \At\ reaches 1.3% at its 
maximum, however in a region, where a < 0.3 fb, the asymmetry At thus cannot 
be measured. In the |//|-M 2 plane for ip Ml = and ip^ = O.bn we found \ Ae\ < 0.7%. 



■ Production of x 2 X3 

For the process e + e _ — > X2X3 we discuss the decay X3 — * ^X? or the heavier neu- 
tralino which has a larger kinematically allowed region than that for x® ~~ * %Xi- 
Similar to x1 X3 production and decay, the cross section ap(e + e~ — > X2X3) reaches 
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Figure 2.22: Contour lines of o = a P (e + e - -> xlxl) x BR(x° -> ^X?) x BR(Z -> 
(2.22a), and the asymmetry ^ (2.22b) in the \fx\-M 2 plane for <p Ml = 0-5tt, cp^ = 0, 
tan/3 = 10, m = 300 GeV, y/s = 800 GeV and (P e -,P e +) = (0.8, -0.6). The area 
A (B) is kinematically forbidden by m^o + m^o > y^s (m z + m^a > m^o). The gray 
area is excluded by m~± < 104 GeV. 

J A.1 



values up to 50 fb for a beam polarization of (P e -,P e +) = (0.8, —0.6) and that for 
the complete process a = a P (e + e~ — > X2X3) x BR(x° — >■ ^X?) x BR(Z — > ££) attains 
values up to 5 fb in the investigated regions of the |//|-M 2 plane in Fig. 2.22a. 



The asymmetry Ae, Fig. 2.22b, is somewhat larger than that for Xi X3 production 
and decay, and reaches at its maximum 2%. However, it will be difficult to measure 
Ae, since e.g. for = 380 GeV, M 2 = 560 GeV and (fMu^fi) = (0.57r,0), we 
found Se ~ 1, for £ = 500 fb~ . For the hadronic decays of the Z boson we have 
Sb(c) ~ 8(5) and thus Ab( c ) is accessible for xl X3 production. For ip^ = 0.5tt and 
(p Ml =0 we found that \Ai\ < 1% in regions of the \fi\-M 2 plane where a < 0.5 fb, 
and \Ae\ < 0.4% in regions where a < 5 fb. 
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2.3.5 Summary of Section 2.3 

We have analyzed CP sensitive observables in neutralino production e + e~ — »■ x°i X°j 
and the subsequent two-body decay of one neutralino into a Z boson xl — > x° Z, 
followed by the decay Z — > for £ = e, /x, r, or Z — > gg with g = c, 6. The 
CP sensitive observables are defined by the vector component V2 of the Z boson 
density matrix and the CP asymmetry Ae( q ), which involves the triple product 
%{q) — Pe- • (Pe(q) x Pi(q))- The tree level contributions to these observables are 
due to correlations of the neutralino x° spin and the Z boson spin. In a numer- 
ical study of the MSSM parameter space with complex Mi and \x for X2X2/ 
X1X3 an d X2X3 production, we have shown that the asymmetry Ae can go up to 
3%. For the hadronic decays of the Z boson, larger asymmetries are obtained with 
A c(b) ~ 6.3(4.5) x A t . 



Chapter 3 



CP violation in production and decay 

of charginos 



Overview 



We study chargino production with longitudinally polarized beams e + e~ — > xt Xj 
with the subsequent leptonic decay of one chargino xt ~^ ?< + vt for i = e, /x, r [48]. 
This decay mode allows the definition of a CP asymmetry which is sensitive to 
the phase (p^ and probes CP violation in the chargino production process. For 
chargino decay into a W boson xt —> W + Xn [49], CP observables can be obtained 
which are also sensitive to tp Ml . We present numerical results for the asymmetries, 
W polarizations, cross sections and branching ratios at a linear electron-positron 
collider with yfe = 800 GeV. 
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3.1 CP asymmetry in chargino production and decay 
into a sneutrino 

We study chargino production 

e + + e- -> xt + Xj; i,j = 1,2, (3.1) 

with longitudinally polarized beams and the subsequent two-body decay of one 
of the charginos into a sneutrino 

Xf^e + + u e ; € = e,Ai,r. (3.2) 
We define the triple product 

% = (Pe- x P x +) • Pe (3.3) 
and the T-odd asymmetry 



T a(T, > 0) - <x(T, < 0) 
At a{% > 0) + cr(7^ < 0) ' 



(3.4) 



of the cross section a for chargino production (3.1) and decay (3.2). The asymmetry 
Aj is not only sensitive to the phase y? M , but also to absorptive contributions, which 
are eliminated in the CP asymmetry 

A e = l(Aj-Aj), (3.5) 

where Aj is the CP conjugated asymmetry for the process e + e~ — > XiXj\ Xi ~^ 
l~ vi. In this context it is interesting to note that in chargino production it is not 
possible to construct a triple product and a corresponding asymmetry by using 
transversely polarized e + and e beams [17,50], therefore, one has to rely on the 
transverse polarization of the produced chargino. 
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3.1.1 Cross section 



For the calculation of the cross section for the combined process of chargino pro- 
duction (3.1) and the subsequent two-body decay of xt (3-2), we use the spin- 
density matrix formalism as in [27,51]. The amplitude squared, 

|T| 2 = |A(x+)| 2 £ Pp(xt) XA < Pnixth^ 0-6) 

is composed of the (unnormalized) spin-density production matrix pp{xt)r de- 
fined in (D.8), and the decay matrix PD{xt)> defined in (D.36), with the helicity 
indices A, and A • of the chargino. Inserting the density matrices into (3.6) leads to 

|T| 2 = 4|A(x+)| 2 (P£ + £^), (3.7) 

where we sum over a. The cross section and distributions are then obtained by 
integrating |T| 2 over the Lorentz invariant phase space element rfLips, defined 
in(B.22): 

1 

da = — |T| 2 rfLips(s;p x -,^,p^). (3.8) 



3.1.2 CP asymmetries 

Inserting the cross section (3.8) into the definition of the asymmetry (3.4) we obtain 

T J Sign[T,] |T| 2 rfLips J Sign^] Ejj]| rfLips 

1 J|T| 2 ciLips JPDdLips ' { ' 

where we have already used the narrow width approximation for the chargino 
propagator. In the numerator of (3.9) only the CP sensitive contribution EpE|, 
from chargino polarization perpendicular to the production plane remains, since 
only this term contains the triple product % = (p e x p +) • p e (3.3). In the denom- 
inator only the term PD remains, since all spin correlations SpEJ, vanish due to 
the integration over the complete phase space. 

The coefficient Ep is non-zero only for production of an unequal pair of charginos, 
e + e — > xfx2> an d obtains contributions from Z-exchange and Z-v interference 
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only, see (D.17). The contribution to T? P from Z-exchange, see (D.27), is non-zero 
only for tp^ ^ 0, it, whereas the Z-v interference term, see (D.28), obtains also ab- 
sorptive contributions due to the finite Z-width which do not signal CP violation. 
These, however, will be eliminated in the asymmetry At (3.5). 

For chargino decay into a tau sneutrino, xt ~^ t + v t , the asymmetry Aj oc ( | V^i | 2 — 
l^i2| 2 )/(|Vii| 2 + Y?\U i2 \ 2 ) is reduced, which follows from the expressions for D 
and Y? D , given in (D.39) and (D.40). 



3.1.3 Numerical results 

We present numerical results for the asymmetries Ae (3.5), for I = e, \x, and the 
cross sections o = o>(e + e~ — > X1X2 ) x BR(x| — > £ + vi)- We study the depen- 
dence of the asymmetries and cross sections on the MSSM parameters fj, — e lipM , 
M 2 and tan (3. We choose a center of mass energy of y/s = 800 GeV and longitu- 
dinally polarized beams with beam polarizations (P e -,P e +) = (—0.8, +0.6), which 
enhance v e exchange in the production process. This results in larger cross sections 
and asymmetries. 

We study the decays of the lighter chargino xt- For the calculation of the chargino 
widths T x + and the branching ratios BR(xt — >■ i^vi) we include the following 
two-body decays, 

Xt -> W + X° n , e + U e , l^Vp,, T + U T , l\v e , [l\v^ f^r, (3.10) 

and neglect three-body decays. In order to reduce the number of parameters, we 
assume the relation \Mi\ = 5/3 M 2 tan 2 $w For all scenarios we fix the sneutrino 
and slepton masses, m^, = 185 GeV, £ = e,fi, r, = 200 GeV, i = e,fi. These 
values are obtained from the renormalization group equations (A.57) and (A.58), 
for M 2 = 200 GeV, m = 80 GeV and tan (3 = 5. In the stau sector, see Appendix 
A.3.3, we fix the trilinear scalar coupling parameter to A T = 250 GeV. The stau 
masses are fixed to m fl = 129 GeV and m f , 2 = 202 GeV. 

In Fig. 3.1a we show the contour lines of the cross section for chargino production 
and decay a = ap(e + e~ — > X1X2) x BR(xi" — > £ + i>e) in the M 2 -(p^ plane for = 
400 GeV and tan/5 = 5. The production cross section o>(e + e~ — > X1X2 ) can attain 
values from 10 fb to 150 fb and BR(xf — > t^vi), summed over t — e, /j,, can be 
as large as 50%. Note that o is very sensitive to y? M , which has been exploited in 
[17, 18] to constrain cos(<y? M ). 
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Figure 3.1: Contour lines of cr = o>(e + e — > X1X2 ) x BR(xi~ 



9+r, 



Vi), summed 



over £ = e, fi, (3.1a), and the asymmetry At for £ = e or \i (3.1b), in the M 2 - 
<Pn plane for \fi\ = 400 GeV, tan/3 = 5, = 185 GeV, y/s = 800 GeV and 
(P e -,P e +) = (—0.8,0.6). The gray area is excluded by m x ± < 104 GeV. The area 
A is kinematically forbidden by mp t + m x o > m x +- The area B is kinematically 
forbidden by + m v - > \fs. 



The M 2 -(ff 1 dependence of the CP asymmetry Ae for £ = e or /x is shown in Fig. 3.1b. 
The asymmetry can be as large as 10% and it does, however, not attain maximal 
values for ip^ = ±0.57r. The reason is that At is proportional to a product of a CP- 
odd (Sp) and a CP-even factor (El,), see (3.9). The CP-odd (CP-even) factor has 
as sine-like (cosine-like) dependence on ip^. Thus the maximum of At is shifted 
towards cp^ = ±tc in Fig. 3.1b. Phases close to the CP conserving points, p^ = 
0, ±7r, are favored by the experimental upper limits on the EDMs, as discussed in 
Section 1.2. 

For M 2 = 200 GeV, we show the tan dependence of a and At in Figs. 3.2a,b. 
The asymmetry can reach values up to 30% and shows a strong tan j3 dependence 
and decreases with increasing tan/3. The feasibility of measuring the asymme- 
try depends also on the cross section a = a P (e + e~ — > X1X2) x BR(x^ — * £ + ve), 
Fig. 3.2a, which attains values up to 20 fb. 
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Figure 3.2: Contour lines of cr = 0>(e + e — > X1X2 ) x BR(xi~ — > £ + ve), summed 
over £ = e, fi, (3.2a), and the asymmetry Ae for I = e or // (3.2b), in the tan (3- 
cp^ plane for M 2 = 200 GeV, = 400 GeV, m„ e = 185 GeV, ^ = 800 GeV and 
(P e -,P e +) = (—0.8,0.6). The area A is kinematically forbidden by m^+m x o > rn x +. 



For the phase ip^ = 0.9tt and tan/3 = 5, we study the beam polarization depen- 
dence of Ae, which can be strong as shown in Fig. 3.3a. An electron beam polariza- 
tion P e - > and a positron beam polarization P e + < enhance the channels with 
v e exchange in the chargino production process. For e.g. (P e -,P e +) = (—0.8, 0.6) 
the asymmetry can attain —7%, Fig. 3.3a, with <j P (e + e~ — > X1X2) ~ 10 fb and 
BR(xt — > £ + v t ) ~ 50%, summed over £ = e,fi. The cross section o = a P (e + e~ — > 
X1X2) x BR(x| -> ^) ranges between 2.3 fb for (P e -,P e +) = (0, 0) and 6. 8 fb for 
(P e -,P e +) = (—1,1). The statistical significance of Ae, given by Si = \Ai\y2C ■ a, 
is shown in Fig. 3.3b for £ = 500 fb -1 . We have S<»5 for (P e -,P e +) = (-0.8, 0.6), 
and thus Ae could be accessible at a linear collider, even for y? M = 0.97T, by using 
polarized beams. 
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Figure 3.3: Contour lines of the asymmetry At for £ = e or ji (3.3a), and the signif- 
icance St (3.3b), for e + e~ — > XiX.2'1 Xi ~> ^ + ^e in the P e --P e + plane for y? M = 0.97T, 
I At I = 400 GeV, M 2 = 200 GeV, tan/3 = 5, m 0t = 185 GeV, = 800 GeV and 
£ = 500 fb" 1 . 



3.1.4 Summary of Section 3.1 



We have studied CP violation in chargino production with longitudinally polar- 
ized beams, e + e~ — > xtxj> an d subsequent two-body decay of one chargino into 
the sneutrino xt —> ^ + ^e- We have defined the T-odd asymmetries Aj of the triple 
product (p e - x p-+) • pt- The CP-odd asymmetries At = h{Aj — AJ), where Aj 
denote the CP conjugated of Aj, are sensitive to the phase tp^ of the Higgsino mass 
parameter ji. At tree level, the asymmetries have large CP sensitive contributions 
from spin-correlation effects in the production of an unequal pair of charginos. In 
a numerical discussion for e + e~ — >■ X1X2 production, we have found that At for 
I = e or \i can attain values up to 30%. By analyzing the statistical errors, we 
have shown that, even for e.g. ip^ ~ 0.97T, the asymmetries could be accessible in 
future e + e" collider experiments in the 800 GeV range with high luminosity and 
longitudinally polarized beams. 
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3.2 CP violation in chargino production and decay into 
a W boson 




Figure 3.4: Schematic picture of the chargino production and decay process. 

CP violation in the chargino sector can be studied also by the two-body decay of 
the chargino into a W boson. In contrast to the decay into a sneutrino, the spin 
correlations of chargino and W lead to CP observables, which are also sensitive to 
the phase of Mi. Since these observables have in addition contributions from the 
decay they do not necessarily vanish for e + e~ — > XiXi production. In this mode 
CP violation could be established, even if the production of the heavier pair X1X2 
is not yet accessible. 

We study chargino production 

e + + e" - xt + Xj; i,i = l,2, (3.11) 
with longitudinally polarized beams and the subsequent two-body decay 

xt^W + + x° n . (3.12) 
We define the triple product 



7} = p e ■ (p + x p w ) 



(3.13) 
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and the T-odd asymmetry 

T _ tr(7} > 0) - crjTr < 0) 
^ " a(T I >0) + a(T I <0Y ^ 

with cr the cross section of chargino production (3.11) and decay (3.12). The asym- 
metry A J is sensitive to the CP violating phase </? M . 

In order to probe also the phase (p Ml , which enters in the chargino decay pro- 
cess (3.12), we consider the subsequent hadronic decay of the W boson 

W + ^c + s. (3.15) 

The correlations between the xt polarization and the W boson polarization lead 
to CP sensitive elements of the W boson density matrix. The triple product 

%i = Pe- • (Pc x p s ), (3.16) 

which includes the momenta of the W decay products, and thus probes the W 
polarization, defines a second T-odd asymmetry 



, T a{?n > 0) ~ opn < 0) 
11 a(T II >0) + a(T II <oy 



(3.17) 



Here, a is the cross section of production (3.11) and decay of the chargino (3.12) 
followed by that of the W boson (3.15). Owing to the spin correlations, Aj T has CP 
sensitive contributions from (p^ due to the chargino production process (3.11) and 
contributions due to y? M and (fM 1 from the chargino decay process (3.12). 

The T-odd asymmetries Aj and AJ T have also absorptive contributions from s- 
channel resonances or final-state interactions, which are eliminated in the CP-odd 
asymmetries 

A I = 1 -(Aj-Aj), A II = ±(A] I -A] I ), (3.18) 

where Aj n are the CP conjugated asymmetries for the processes e + e" — > XiXj'i ' 
W'Xn and e+e" -> XiXj\ Xl -»• W~xl\ W~ -»• cs, respectively. 
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3.2.1 Spin density matrix of the W boson 

For the calculation of the amplitudes squared for the combined process of chargino 
production (3.11) and the subsequent two-body decays (3.12) and (3.15) of xf we 
use the same spin-density matrix formalism as in [27,51]. The (unnormalized) 
spin-density matrix of the W boson 

p P (W + ) x ^ = \A(xt)\ 2 £ Pp(xt) XlK PnAxtfjt ( 3 - 19 ) 

is composed of the chargino propagator A(xf), the spin-density production matrix 
Pp(xt)' defined in (D.8), and the decay matrix pDi(xt)r defined in (D.42). The 
amplitude squared for the complete process e + e~ — > xtxj'f xt ~^ W + Xn'> W + — > 
/'/ can now be written 

|T| 2 = |A(^+)| 2 £ Pp(W + )^ PD2 (W + ) KXk , (3.20) 

Afc ,A fc 

with the decay matrix for W decay pd 2 {W + ), defined in (D.45). Inserting the den- 
sity matrices pp{xt) (D-8) and PdAxI) (D.62) into (3.19) gives 

p P {W + ) XkX * = 4 |A(x+)| 2 [(PD 1 + E P E^) 5 XkX * + (P C D 1 + T? P e E a Dl ) (J c ) AfcA * 

+(P cd D l + E P cd E^) {J cd ) Xk % (3.21) 

summed over a, c, d. We have thus decomposed the W production matrix pp(W + ) 
into contributions of scalar (first term), vector (second term), and tensor parts 
(third term). Inserting then p P (W + ) (3.21) and Pd 2 (W + ) (D.63) into (3.20) gives 

|T| 2 = 4|A(x+)| 2 |A(^+)| 2 {3(PP 1 + E P E r jP 2 + 2(P^D 1 + E P c E r J c P 2 

+4[(P cd D 1 + E P cd Z a Di ) cd D 2 - |(P CC D 1 + E P CC E^) dd D 2 ]}. (3.22) 

3.2.2 W boson polarization 

The mean polarization of the W bosons in the laboratory system is given by the 
3x3 density matrix < p(W + ) >, obtained by integrating (3.21) over the Lorentz 
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invariant phase-space element dLips(s;p x -,p x o,p w ), see (B.39), and normalizing 
by the trace 

< p(W + ) x ^ >= S Pp(W + ) x ^ dLij>s = l5 x k y y (JC) A fc A' fc+T {jcd) x k x> k (3<23) 

The components V c of the vector polarization and T cd of the tensor polarization are 

_ J(P c D 1 + £^£ a Di )rfLips 
Vc ~ 3/PAdLips ' ( j 

j{P cd D l + ££ cd S^) rfLips 
3 J aLips 

where we have already used the narrow width approximation for the chargino 
propagator. The density matrix elements in the helicity basis (B.38) are given by 

<p{W + )- > = \-V 3 + T 33 , (3.26) 
< p(W + ) 00 > = -2 T 33 , (3.27) 
<p(W + )-°> = ^(V 1 + iV 2 )-V2(T 13 + iT 23 ), (3.28) 

< p(W + )- + > = T n - T 23 + 2tT l2 , (3.29) 
<p(W + ) 0+ > = -±(V 1 +tV 2 ) + V2(T n + iT 23 ), (3.30) 

where we have used Tn + T 22 + T 33 — — ». 



3.2.3 T-odd asymmetries 

From (3.21) we obtain for asymmetry Aj (3.14): 

T _ fSign[T I }Ti{p P (W+) x ^} dLips _ /SignpTj^S^dLips 

7 / Tr {p P ( ) x ^ } dLips /PDi dLips ' K ' ' 

with c?Lips(s;p x -,p x o,j9 W ) given in (B.39), where we have used the narrow width 
approximation. In the numerator of (3.31), only the spin correlations £f> T? Dl per- 
pendicular to the production plane remain, since only this term contains the triple 
product 7} = p e - • (p + x pw)- In the denominator only the term PD 1 remains, and 
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all spin correlations vanish due to the integration over the complete phase space. 
Note that Aj oc Y? Dl oc (|0^| 2 - |0^| 2 ), see (D.59), and thus Aj may be reduced 
for 1 Rj \0^\. Moreover, Aj will be small for m 2 + — m 2 ~ 2m^, see (D.59). 

For the asymmetry Ajj (3.17), we obtain from (3.22): 

T /Sign[T„]|r| 2 rfLips JSign[T JJ ]2S» C D 2 rfLips 

11 J|T| 2 rfLips /3PDiD 2 dLips ' { ' ' 

with dLips(s;p x -,p x o,Pf>,Pf), defined in (B.40), where we have used the narrow 

width approximations for the propagators. In the numerator only the vector part 
of |T| 2 remains because only the vector part contains the triple product T n = p e ■ 
(p c x p 3 ). In the denominator the vector and tensor parts of |T| 2 vanish due to 
phase-space integration. Owing to the correlations E P c T, a Di between the xt an d 
the W boson polarization, there are contributions to the asymmetry Ajj from the 
chargino production process (3.11), and/ or from the chargino decay process (3.12). 
The contribution from the production is given by the term with a = 2 in (3.32) 
and it is proportional to the transverse polarization of the chargino perpendicular 
to the production plane Ep, see Appendix D.1.2. For the production of a pair of 
equal charginos, e + e~ — > xtxi> we have Ep = 0. The contributions from the decay, 
given by terms with a = 1, 3 in (3.32), are proportional to 

C E^ C D 2 D -2g*m xl lm(0^0^)(t c w ■ P/ ) w p» + f w t% , (3.33) 

see last term of (D.60), which contains the e-tensor. Thus Ajj can be enhanced 
(reduced) if the contributions from production and decay have the same (opposite) 
sign. Note that the contributions from the decay vanish for a two-body decay 
of the chargino into a scalar particle instead of a W boson. In order to measure 
Ai the momentum of xt* i- e - the production plane, has to be determined. This 
could be accomplished by measuring the decay of the other chargino xj- F° r the- 
measurement of An, the flavors of the quarks c and s have to be distinguished, 
which will be possible by flavor tagging of the c-quark [30,32]. In principle, for 
the decay W — > u d also an asymmetry similar to An can be considered, if it is 
possible to distinguish between the u and d jet, for instance, by measuring the 
average charge. 

3.2.4 Numerical results 

We study the dependence of Ai, An, and the density matrix < p(W + ) > on the 
MSSM parameters fi = e 1 ^, M 1 = \Mi\e l(pM i, tan/3 and the universal scalar 
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mass parameter m . The feasibility of measuring the asymmetries depends also 
on the cross sections o = a P (e + e~ — > xtxj) x BR(x+ — > W^+x^) x BR(VV + — > cs), 
which we will discuss in our scenarios. We choose a center of mass energy of y/s = 
800 GeV and longitudinally polarized beams with (P e -, P e +) = (—0.8, +0.6). This 
choice enhances sneutrino exchange in the chargino production process, which re- 
sults in larger cross sections and asymmetries. For the calculation of the branching 
ratios BR(x+ — > W + xV) and widths T x +, we include the two-body decays 

xt -> W + x° n , e£v e , Ji\v^ f+ 2 v T , e + z> e , /i + z> M , t + u t , (3.34) 

and neglect three-body decays. For the W boson decay we take the experimental 
value BR(Vy + — > cs) = 0.31 [45]. In order to reduce the number of parameters, 
we assume the relation \M±\ = 5/3 M 2 tan 2 9w and use the renormalization group 
equations for the slepton and sneutrino masses, see Appendix A.3.4. In the stau 
sector, see Appendix A.3.3, we fix the trilinear scalar coupling parameter A T = 
250 GeV. 

■ Production ofxt Xi 

For the production e + e — >■ xtxT °f a P a i r °f charginos the polarization perpen- 
dicular to the production plane vanishes, and thus Ai — 0. However, „4// need 
not to be zero and is sensitive to and (p Ml , because this asymmetry has con- 
tributions from the chargino decay process. For (</?Mi, <Ph) = (0.5n, 0) we show in 
Fig. 3.5a the \/jl\-M 2 dependence of An, which can reach values of 5% — 7% for 
M 2 > 400 GeV. We also studied the dependence of An in the \fi\-M 2 plane. For 
V?Mi = 0, y? M = 0.l7r(0.57r) and the other parameters as given in the caption of 
Fig. 3.5, we find \A n \ < 2%(7%). 

In Fig. 3.5b we show the contour lines of the cross section o = a P (e + e~ — > XiXi ) x 
BR(x^ -> W^ + X?) x BR(Vy+ -> cs) in the |//|-M 2 plane for ((p Ml ,<Pp) = (0.5tt,0). 
The production cross section <rp(e + e _ — > x^x^) reaches more than 400 fb. For our 
choice of m = 300 GeV, xt — ^ is the only allowed two-body decay channel. 

In Fig. 3.6a we plot the contour lines of An for M = 350 GeV and M 2 = 400 GeV 
in the ip^-ip^ plane. Fig. 3.6a shows that An is essentially depending on the sum 
V/i+^Mi- However, maximal phases of y? M i = ±0.57r and v?^ = ±0.57rdonotleadto 
the highest values of \ An\ > 6%, which are reached for ((p Ml ^n) ~ (±0.87r, ±0.67r). 
The reason for this is that the spin-correlation terms Sp c SJ, i C D 2 in the numerator 
of An (3.32) are products of CP-odd and CP-even factors. The CP-odd (CP-even) 
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Figure 3.5: Contour lines of the asymmetry An (3.5a) and cr = tjp(e + e — > XiXi ) x 
BR(x^ -> W + x\) x BR(iy+ -> cs) (3.5b), in the |/i|-M 2 plane for {<p Ml ,<pJ = 
(0.5?r, 0), tan/3 = 5, m = 300 GeV, y/s = 800 GeV and (P e -,P e+ ) = (-0.8, 0.6). The 
area A (B) is kinematically forbidden by m x + + m x - > y/s (m w + m x o > m x +). The 
gray area is excluded by m x ± < 104 GeV. 



factors have a sine-like (cosine-like) phase dependence. Therefore, the maximum 
of the CP asymmetry An may be shifted to smaller or larger values of the phases. 
In the (p/jr-ipMi region shown in Fig. 3.6a the cross section a = a P (e + e~ — > XiXi) x 
BR(x^ -> W+xD x BR(W + -> cs) with BR(^~ -> = 1 does not depend on 

(pMx an d ranges between 74 fb for = and 66 fb for (p^ = tt. 

In Fig. 3.6b we show the contour lines of the significance Sn = \An\V2£ ■ a, de- 
fined in (1.18). For C = 500 fb -1 and for e.g. (y?A/ i; </v) ~ ( n i 0.l7r) we have 5*// « 8 
and thus Azr should be measured even for small ip^. 
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Figure 3.6: Contour lines of and the asymmetry An (3.6a) and the significance 
S n (3.6b) for e+e~ -> XiXi', Xi -»• W + x°i, W + -> cs, in the p lx -p Ml plane for 
|/x| = 350 GeV, M 2 = 400 GeV, tan (5 = 5, m = 300 GeV, = 800 GeV, (P e - , P e+ ) = 
(—0.8, 0.6) and C = 500 fb _1 . In the gray shaded area of Fig. 3.6b we have Sn < 5. 



In Figs. 3.7a,b we show the tan/3-m dependence of An and o for (tp Ml ,p fJi ) = 
(0.77T, 0). The asymmetry is rather insensitive to m and shows strong depen- 
dence on tan (3 and decreases with increasing tan (3 > 2. The production cross sec- 
tion o>(e + e~ — > XiXi) increases with increasing m and decreasing tan/3. For 
m < 200 GeV, the branching ratio BR(xf — »■ < 1/ since the decay channels 

of Xi" into sleptons and/ or sneutrinos open, see (3.34). 

In Fig. 3.8a we show the tp^ dependence of the vector (Vi) and tensor (T^) com- 
ponents of the density matrix < p(W + ) > for (p Ml = it. In Fig. 3.8b we plot their 
V?Mi dependence for p^ = 0. In both figures, the element V 2 is CP-odd, while 
T\z, Tu, T 2 2 and Vi, V 3 show a CP-even behavior. As discussed in Appendix D.3, 
the components T n and T 22 are almost equal and have the same order of magni- 
tude as Vi and V 3 , whereas T 12 , |T 23 | < 10~ 5 are small. For CP conserving phases 
(p Ml ,Pfi) = (0, 0) the density matrix reads 



< p(W + ) > 



(<P- 
<p G 



> < P 
> 



< P °° > 



> < P 



-+ 



> 



> < p +0 > < p ++ > 



( 0.200 -0.010 -0.001 \ 



-0.010 
-0.001 



0.487 
0.137 



0.137 
0.313 j 



, (3.35) 
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Figure 3.7: Contour lines of the asymmetry An (3.7a) and a = o>(e + e — > XiXi ) x 

BR(x^ -> x BR(W + -> cs) (3.7b), in the tan/5-m plane for ((p Ml ,<Pp) = 

(0.7tt, 0), M 2 = 400 GeV, = 350 GeV, ^ = 800 GeV and (P e -,P e+ ) = (-0.8, 0.6). 



for M 2 = 400 GeV, |//| = 350 GeV, tan/5 = 5, m = 300 GeV, ^ = 800 GeV, 
(P e -,P e +) = (—0.8, 0.6) For CP violating phases, e.g. (y9 Ml , y? M ) = (0.7 it, 0) and 
the other parameters as above, the density matrix has imaginary parts due to a 
non-vanishing V 2 : 



0.219 -0.010 + 0.025i 0.002 
<p(W + ) >= | -0.010 - 0.025* 0.405 0.171 + 0.025i 

0.002 0.171 -0.025* 0.376 



(3.36) 



Imaginary parts of the density matrix are thus an indication of CP violation. 



■ Production of xt X2 



For the production of an unequal pair of charginos, e + e — > X1X2 > their polar- 
ization perpendicular to the production plane is sensitive to the phase y? M , which 
leads to a non-vanishing asymmetry Ai (3.31). We will study the decay of the 
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Figure 3.8: Dependence of vector (V*) and tensor (T^ ) components of the den- 
sity matrix on p^ (3.8a) and on (3.8b), for e + e~ — > x^xr; Xi" ~~ *■ W + X?/ ror 
= 350 GeV, M 2 = 400 GeV, tan/3 = 5, m = 300 GeV, y^ = 800 GeV and 
(P e -,P e+ ) = (-0.8,0.6). 



lighter chargino Xi" - ► W + Xi- F° r I ^2 1 = 250 GeV and y9 Ml = 0, we show in 
Fig. 3.9a the \fi\-ipn dependence of Ai, which attains values up to 4%. Note that 
Ai is not maximal for ip^ = 0.5n, but is rather sensitive for phases in the regions 
p^ e [0.77r,7r] and p^ G [— 0.77T, — tt]. As discussed in Section 1.2, values of y? M 
close to the CP conserving points p^ = 0, ±ir are suggested by EDM analyses. For 
Pf, = 0.97T and = 350 GeV the statistical significance is Si = \Ai\V2L - a pz 1.5 
with C = 500 fb _1 . Thus Ai could be measured at a confidence level larger than 
68% (Si = 1). 

In Fig. 3.9b we show contour lines of the corresponding cross section a = a P (e + e~ - 
X1X2 ) x BR(x^ — ► W + X?) i 11 t ne l/^hw plane for the parameters as above. The 
cross section shows a CP even behavior, which has been used in [17, 18, 50] to con- 
strain cos (pp. In our scenario we have considered the decay of the lighter chargino 
Xi — > W + Xi since for our choice m = 300 GeV we have BR(x^ — > W /+ Xi) = 1- 
For the decay of X2" / one would have to take into account also the decays into the Z 
boson and the lightest neutral Higgs, which would reduce BR(x^ — ► W /+ Xi) ~ 0.2. 
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Figure 3.9: Contour lines of the asymmetry Ai (3.9a) and a = a P (e + e — > x^x 2 ) x 
BR(x^ -> VT+x?) (3.9b), in the |//|-^ plane for y> Ml = 0, M 2 = 250 GeV, tan/5 = 5, 
m = 300 GeV, js = 800 GeV and (P e -,P e +) = (-0.8,0.6). The area A (B) is 
kinematically forbidden by m x + + m x - > \fs (mw + m x a > m x +)- The gray area is 
excluded by m x ± < 104 GeV. 



The asymmetry An is also sensitive to the phase ipMi- We show the f^-^Mi de- 
pendence of An, choosing the parameters as above, in Fig. 3.10a. In Fig. 3.10b we 
show the contour lines of the significance Sn = \An\y2C ■ a for £ = 500 fb _1 . For 
(ipMntPn) ~ I 71 "? 0.l7r) we have Sn « 2.4 and thus An could be accessible even for 
small phases by using polarized beams. 



3.2.5 Summary of Section 3.2 

We have analyzed CP sensitive observables in chargino production, e + e~ — > xfxj, 
with subsequent two-body decay, xt —> W + Xn- We have defined the CP asymme- 
try Ai of the triple product p e • (p^+ x pw)- In the MSSM with complex parameters 
/i and Mi, we have shown that Ai can reach 4% and that even for cp^ « 0.97T the 
asymmetry could be accessible in the process e + e _ — > ^1^X2 • Further we have an- 
alyzed the CP sensitive density-matrix elements of the W boson. The phase <p Ml 
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Figure 3.10: Contour lines of and the asymmetry An (3.10a) and the significance 
Sn (3.10b) for e + e~ — > XiX.2'1 Xi ~^ W + Xi! W + — »■ cs, in the ip^-ipMi plane for 
= 350 GeV, M 2 = 250 GeV, tan (3 = 5, m = 300 GeV, Vs = 800 GeV, (P e - ,P e +) = 
(-0.8, 0.6) and £ = 500 fb _1 . In the gray shaded area of Fig. 3.10b we have S n < 1. 



enters in the decay xt ~^ W + Xn due to correlations of the chargino and the W 
boson spins, which can be probed via the hadronic decay W + — > cs. Moreover the 
triple product p e • (p c x p^) defines the CP asymmetry An, which can be as large 
as 7% for XiXi or X1X2 production. By analyzing the statistical errors of Ai and 
An we found that the phases y? M and </?Mi could be strongly constrained at a future 
e + e~ collider with y/s = 800 GeV, high luminosity and longitudinally polarized 
beams. 



Chapter 4 



CP violation in sf ermion decays 



For the sfermion decays 

+ X°->£ + X°l Z->£ + I £ = e,n,r, (4.1) 
+ + Z^q + q, q = b,c, (4.2) 

schematically shown in Fig. (4.1), the triple product of the momenta of the outgo- 
ing leptons 

% = Pf ■ (Pi x Pi), (43) 
and that of the outgoing quarks, 

T q = Pf {P q x Pf), (4-4) 




Figure 4.1: Schematic picture of the sfermion decay process. 
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define T-odd asymmetries 

T r(T^>o)-r(T £ , g <o) 
^ r(T £ , g >o) + r(T, i(? <o) ^ 

of the partial sfermion decay width Y for the process (4.1). The asymmetries [52] 
are sensitive to correlations between the x°j polarization and the Z boson polariza- 
tion, which are encoded in the momenta of the final leptons or quarks. The corre- 
lations thus would vanish if a scalar particle in place of the Z boson is exchanged. 
The tree-level contribution to the asymmetries (4.5) are proportional to the imag- 
inary part of a product of the x'j'Z'Xi couplings. However, they are not sensitive 
to the phase (pA f of the trilinear scalar coupling parameter Af, since the decay 
/ — * fx°j (4-1) is a two-body decay of a scalar particle. As an observable in the 
process (4.1) which is sensitive to tpA f , one would have to measure the transverse 
polarization of the fermion /, which is possible for f = r,t [53]. Also three-body 
decays of the sfermion / can be studied alternatively [54]. 

In the numerical study we estimate the event rates necessary to measure the asym- 
metries, which can reach 3% for leptonic decays t x — »■ r xi ? an d 20% for semi- 
leptonic decays like t x — > r x\ b b. 

The triple product (4.3) was proposed in [55] and the size of the asymmetry was 
calculated for the decay \x — > \i x\ — > Xi f 1 ? ^> however, for a specific final state 
configuration only. We extend the work of [55] by calculating the asymmetries (4.5) 
in the entire phase space. 



4.1 Sfermion decay width 

For the calculation of the amplitude squared of the subsequent two-body decays 
of the sfermion (4.1), we use the spin-density matrix formalism of [27]: 

m 2 =|A(^°)| 2 |A(zT)| 2 £ f PoAfXx'Pn^ZiPosi 2 )^- (4-6) 

The amplitude squared is composed of the propagators A(x°), A(Z), the unnor- 
malized spin density matrices po 1 (/), Pd 2 (x°) and Pd 3 (Z), with the helicity indices 
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Xj, X'j of the neutralino and/or the helicity indices \' k of the Z boson. Intro- 
ducing a set of spin four- vectors s a , a = 1,2,3, for the neutralino x°i> see Ap- 

Xj 

pendix (B.49), the density matrices can be expanded in terms of the Pauli matrices 



"3" j ' 



b fiu 



(4.7) 
(4.8) 
(4.9) 



The polarization vectors e^ k of the Z boson obey p^e^ = and the completeness 
relation J2x k ^ k *^l k = —g^u+Pz^Pzv/ , m 2 z- The expansion coefficients of the density 
matrices (4.7)-(4.9) are 



D 1 



D 



2 pa 



H\ 2 + U 3 \ 2 )(P f i-P x o), 

±m x o(\al j \ 2 -\b{ j \ 2 )(p f/ .s%), 



{g pa [2 Re(0'(fO'lf* )m x? m x o - (|0^| 2 + \0»*\ 2 )(p# ■ p x o) 



cos 2 9w 

+ (\0'^\ 2 + \0'[f | 2 )(p x> p x?CT +P x o /x?p )}, 



cos 2 9 W 



{2 m x olm(0»?0»f)(p x o p s\ - p x> «V) 



-w^o^(i^i 2 + iorfn 

+ 2 W p^^m x oRe(0 / 1 'j : / 1 / f )}, 



2 5 s 



cos 2 # 



w 



{ - <T{L) + R})(p f ■ pj) + (p p f pj + P P jP° f ){L) + R}) 



i{R) 



L)) e pupv 



(4.10) 
(4.11) 



(4.12) 



(4.13) 



(4.14) 



with e i23 = 1 and the couplings as defined in Appendix A.2. The negative sign 
in (4.11) holds for the decay of a negatively charged sfermion. In (4.10) and (4.11), 
/' denotes the fermion from the first decay /—>/'%!■ in (4.1). Inserting the density 
matrices (4.7)-(4.9) in (4.6), we obtain for the amplitude squared 

\T\ 2 = 2 |A(x?)| 2 |A(zT)| 2 {L>i D 2pa + ^ a D , 2pa }D p °. (4.15) 
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The / decay width for the decay chain (4.1) is then given by 

r(/~ ^f'Xiff) = ^:J \T\ 2 dUps(mf,p r ,p x o,p f ,p f ), (4.16) 
with the phase-space element ciLips defined in Appendix B.3.2. 



4.2 T-odd asymmetry 

In the following we present in some detail the calculation of the T-odd asymme- 
try (4.5) for the slepton decays £ — > £ x°j — > ( Xi Z ~^ ? Xi f /• The replacements 
that must be made to obtain the asymmetry for q decays are obvious. From (4.5) 
and (4.15) we find 



where we have already used the narrow width approximation for the propagators. 
In the numerator we have used / Sign[7^ ;? ]Di Z^po- -DjT^Lips = and in the denom- 
inator / E a D2pa D$ a dUps = 0. Among the spin correlation terms Y? Dl E a D2pa D^ 
only those contribute to Aj q , which are proportional to the triple product T^ q 

^d 2P X => 32 m x M{0'[]0'[f){R) - L))e^ p x> P/ „ P/ „, (4.18) 

see first term of (4.13) and the last term of (4.14). From the explicit representations 
of the neutralino spin vector (B.49) and the lepton momentum vector (B.46), we 
find in the sfermion rest frame (pe ■ s a ) = for a = 1, 2 so that Eg = in (4.11). 

Thus only E|> 2/0(T contributes and the momentum dependent part of (4.18) can be 
written as 



e pa ^ p x o p s x o a p f „ p fv = m x o p e • ( P/ x P/ ), (4.19) 
with p = p/|p|. 



4.3 Numerical results 



75 



The dependence of Aj on the t\ct-Xj couplings on the Z-f-f couplings 

y/L 



L f , R f and on the Z-x\-x] couplings 0"f ' R follows from (4.17) and (4.18): 



\ n t 12 _ \U 12 7-2 _ d2 

A f * ffeM Im( °" °« '>• (420) 

K/ + Wkj? L f + K f 

lj 12 I fo £ |2 

Due to the first factor f f the asymmetry will be strongly suppressed for 

\ a kj\ +l tl fcjl 

\ a kj\ ~ l^ijl an< i maximally enhanced for vanishing mixing in the slepton sector. 
Due to the second factor L 2 +fi 2 , the asymmetry Aj^ for hadronic decays of the Z 
boson is larger than the asymmetry A] for leptonic decays: 



■4c) = L A + % , IT A J * 6.3 (4.5) x Aj. (4.21) 
L i ~ K i L b(c) + K b(c) 

For the measurement of Aj the charges and the flavors of / and / have to be dis- 
tinguished. For / = e, pi this will be possible on an event by event basis. For f — t 
one has to take into account corrections due to the reconstruction of the r momen- 
tum. For f — b,c the distinction of the quark flavors should be possible by flavor 
tagging [30,31]. However, in this case the quark charges will be distinguished 
statistically for a given event sample only. 



4.3 Numerical results 



We assume that t x is the lightest sfermion and study the decay chain f\ — > rx°; x°j ~ 
X°\Z\ Z — > 11 for the two cases fi — > t\\ and f\ — > separately. We present nu- 
merical results for the T-odd asymmetry „4j (4.5) and the branching ratios BR(ri — > 
rx!«):= BR(fi -> r x°) x BR(x° -> x BR(Z -> The size of the asym- 
metry Aj c for hadronic decays may be obtained from (4.21). 

The relevant MSSM parameters are y? M , \Mi\,tp Ml ,M 2 , tan (3, \A T \,^ AT ,m fl and 
m f2 . We fix tan f3 = 10, \A T \ = 1 TeV, <^ At = 0, m fl = 300 GeV, m fa = 800 GeV and 
use the relation \Mi\ = 5/3 M 2 tan 2 6 W in order to reduce the number of parame- 
ters. 
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Table 4.1: Masses and widths for various combinations of <£> M and y>Mi/ for |/x 
300 GeV, M 2 = 280 GeV, tan/3 = 10, A T = 1 TeV, m fl = 300 GeV, m f2 = 800 GeV. 





{ Pm 1 




m x o, 


X s ' 


x o [GeV] 


T x o [MeV] 


I* [MeV] 








135, 


234, 


306, 


358 


4.06 


527 





7T 

2 


137, 


233, 


308, 


357 


1.79 


550 





7T 


138, 


231, 


309, 


356 


0.09 


573 


7T 

2 





137, 


239, 


307, 


353 


5.43 


487 


7T 

2 


7T 

2 


138, 


238, 


309, 


352 


2.89 


511 


7T 

2 


7T 


137, 


237, 


311, 


351 


1.49 


529 


7T 





138, 


245, 


309, 


347 


7.25 


448 


7T 


7T 

2 


137, 


244, 


311, 


346 


5.78 


466 


7T 


7T 


136, 


243, 


313, 


345 


4.32 


484 



For the calculation of the branching ratios BR(fi — > r x°j) and BR(y° — > ^X?)/ we 
include the decays 

ri^rx°, X~^r, (4.22) 
X° - Zx?, xM, m = 1, 2, n < j. (4.23) 

We fix the Higgs mass parameter m A = 800 GeV so that the decays of the neu- 
tralino into charged Higgs bosons —> XmH T , as well as decays into heavy neu- 
tral Higgs bosons x°j —> Xn -^2,3 are forbidden. The decays via sleptons x° — > £ £ 
are forbidden due to our assumption that t x is the lightest sfermion. 



4.3.1 Decay chain via f\ — > r X2 

We study for the decay chain fx — > t\\] x\ —> Xi^; Z —> ££ ior £ = e, //, r. In 
Fig. 4.2a we show the contour lines for the branching ratio BR(ti ^ t Xi £ £) = 
BR(fi -> r x°) x BR(x° -> Zxi) x BR(Z -> £ f), summed over £ = e, /1, r, in 
the M 2 -|//| plane for </? Ml = n/2 and = 0. The t\-t~x\ coupling \d[ 2 \ is larger, 
which implies a larger BR(fi — > r if we choose > M^. We use the usual 
notation M% = M Rf , M L = M Lf , see (A.46) and (A.47). The ordering > M L 
is suggested in some scenarios with non-universal scalar mass parameters at the 
GUT scale [56]. Furthermore, in (A.46) and (A.47) one could have M fRR > M fLL in 
extended models with additional D- terms [57]. 



4.3 Numerical results 
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than for M% < Mg. In a large region of the parameter space BR(x° — ► Zx\) = 1. 
The asymmetry ^4 J is shown in Fig. 4.2b. The dependence of A J on M 2 and is 
dominantly determined by lm.{0'[ 2 0'[ 2 *), as expected from (4.20). 

In Fig. 4.3 we show the </?Mi and y? M dependence of BR(fi — »■ t x\ ££) and of ^ 
for \n\ = 300 GeV and M 2 = 280 GeV. For these parameters, we also give in Table 
4.3.1 the neutralino masses and the total neutralino and stau widths for various 
phase combinations. Note that maximal phases ip^, (p Ml = ±7r/2 do not lead to the 
highest value of Aj, since the asymmetry is proportional to a product of CP even 
(E 3 Dl ) and CP odd terms {T? D2pa D^ ), see (4.17). We give a lower bound on the 
number N of fi's to be produced at a linear collider, in order to measure Aj at la. 
We estimate N = [(Aj) 2 x BR] -1 from the relative statistical error of the asymmetry, 
see (1.16), with BR = BR(ti -^r^f I). For the point y? M = -n/2 and y?Mi = 7r/2, 
marked by • in Fig. 4.3, BR « 2.5% and \Aj\ « 3%, so that iV « 4.4 x 10 5 . For the 
decay h -> 66x?r, however, BR « 3.6% and |^| « 19%, so that only N « 7.7 x 10 2 
fi's are needed. We obtain almost the same results for smaller CP phases y? M = 
and y?Mi = — 0.37T, marked by <g> in Fig. 4.3. In these two examples Aj q should 
be measurable at an e + e _ linear collider with y/s = 800 GeV and an integrated 
luminosity of 500 fb _1 . 



4.3.2 Decay chain via fi — > r x§ 

We discuss the decay chain fi — > rx 3 ; x 3 — > Xi-^j ^ ~^ ££ior £ = e, /i, r. The decay 
fi — > r %3 can be distinguished from fi — > r x° by the r energy. In Fig. 4.4a we show 
the contour lines of BR(n ^tx\££)= BR(fi -> r £°) x BR(x§ -> ^X?) x BR(Z -> 
£ £) in the \fi\-M 2 plane for </? Ml = tt/2 and y?^ = 0. We choose M% < M~ L since the 
fi-r-Xg coupling |a[ 3 | is larger, thus BR(fi — ► r x 3 ) is larger than for > Mg. 
However, the total branching ratio is smaller than for the previous decay chain 
due to the small BR(fi — > r x 3 ) < -75(0.05) in the upper (lower) part of Fig. 4.4a. 

The asymmetry Aj, shown in Fig. 4.4b, vanishes on contours where either \a[ 3 \ = 
\b{ 3 \ or \m(0'l^O'l§*) = 0, see (4.20). Along the contour Aj = in the lower part 
of Fig. 4.4b we have |a[ 3 | = |6[ 3 |, whereas along the contour line in the upper 
part of Fig. 4.4b we have lm(0" 3 O'/jf*) = 0. Furthermore, there is a sign change of 
lm(0'^0'if) between the upper and the lower part of Fig. 4.4b (area A). The first 
factor in (4.20) increases for \fJ,\/M 2 — > 0, since the gaugino component of x 3 gets 
enhanced, resulting in |&[ 3 |/|a[ 3 | — > 0. 

In Fig. 4.5 we show the ifMv W dependence of BR(ri — > r £ £) and *Aj, for = 
150 GeV and M 2 = 450 GeV. Two points of level crossing appear at ((p Ml ,(p ll ) ps 
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(±0.957r, ±0.77r) in Fig. 4.5b. Neutralino masses and the neutralino and stau widths 
are given in Table 4.2. Comparing Fig. 4.3b and Fig. 4.5b, one can see the com- 
mon and strong (p Ml dependence of the asymmetries. In a good approximation 
Sign[„4j] Sign[</? Ml ] in Fig. 4.3b and Sign[.Aj] ps -Sign[y2 Ml ] in Fig. 4.5b, due to 
the different phase dependence of Im(O"|'0i'f *) and lm.(0"^0"f). 



4.4 Summary of Chapter 4 

We have considered a T-odd asymmetry in the sequential decay of a sfermion / — > 
/' X°j —> f Xi Z — > /' Xi //• The asymmetry is sensitive to the phases in the 
neutralino sector. In a numerical study for stau decay fx — > /// we have shown 
that the asymmetry can be of the order of 3% for leptonic final states r x\ an d 
is larger by a factor 6.3 for the semi-leptonic final state r %5 The number of 
produced f's which are necessary to observe the asymmetry is at least of the order 
10 5 for leptonic final states, and 10 3 for semi-leptonic final states, such that the 
phases in the neutralino sector may be accessible at future collider experiments. 
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Figure 4.2: Contour lines of the branching ratio for fi — > y^rU an d asymmetry 
Aj in the /i-M 2 plane for </? Ml = ir/2 and y? M = 0, taking tan/5 = 10, A T = 1 TeV, 
m fl = 300 GeV, m f , 2 = 800 GeV for > M^. The gray areas are kinematically 
forbidden by < m x o + m T (light gray) or m x o < m x o + m z (dark gray). 



Table 4.2: Masses and widths for various combinations of y? M and ipM lf for 
150 GeV, M 2 = 450 GeV, tan/? = 10, A T = 1 TeV, m fl = 300 GeV, m f . 2 = 800 GeV. 









m x o,m x o,m x o [GeV] 


T x o [MeV] 


r fl [MeV] 








128, 


156, 238, 467 


59.0 


362 





7T 

2 


132, 


153, 238, 466 


68.2 


359 





7T 


141, 


145, 238, 466 


75.5 


356 


7T 

2 





131, 


158, 237, 466 


41.5 


356 


7T 

2 


7T 

2 


136, 


154, 237, 466 


49.4 


353 


7T 

2 


7T 


142, 


145, 240, 465 


73.8 


360 


7T 





135, 


159, 236, 465 


27.7 


351 


7T 


7T 

2 


137, 


154, 239, 465 


47.5 


357 


7T 


7T 


143, 


144, 242, 464 


71.0 


364 
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Figure 4.3: Contour lines of the branching ratio for f\ — > xlrii and asymmetry 
Aj in the y? M -y? Ml plane for |//| = 300 GeV, M 2 = 280 GeV, taking tan/3 = 10, 
A r = 1 TeV, m fl = 300 GeV, m f2 = 800 GeV for M^> M~ L . 
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Figure 4.4: Contour lines of the branching ratio for fx — > x5t££ and asymmetry ^4 

in the |//|-M 2 plane for ip Ml = tt/2, y? M = 0, tan /3 = 10, A T = 1 TeV, m fl = 300 GeV, 
rrif 2 = 800 GeV and M% < M^. The area A (B) is kinematically forbidden by 
m^o < m^.o + m z {m fl < m^o + m T ). The gray area is excluded by m x ± < 104 GeV. 
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Figure 4.5: Contour lines of the branching ratio for fi — > x?t£<? and asymmetry Aj 
in the yy-^Mi plane for |//| = 150 GeV, M 2 = 450 GeV, tan/5 = 10, A r = 1 TeV, 
m fl = 300 GeV, m f2 = 800 GeV and M & <M~ L . 



Chapter 5 



Summary and conclusions 



5.1 Summary 

In supersymmetric (SUSY) extensions of the Standard Model (SM), several param- 
eters can be complex. In the neutralino sector of the Minimal extended Supersym- 
metric Standard Model (MSSM), these are the Higgsino mass parameter /j, and the 
gaugino mass parameter M±. In addition, in the sfermion sector also the trilinear 
scalar coupling parameter Af can be complex. The imaginary parts of these pa- 
rameters imply CP violating effects in the production and decay of SUSY particles. 

In this thesis we have analyzed the implications of tp^, <p Ml and tp^ in neutralino 
and chargino production and decay in electron-positron collisions. For the decays 
of sfermions we have analyzed the effects of p^ and tp Mi • We have analyzed T-odd 
and CP-odd asymmetries of triple products of particle momenta or spins. Such 
asymmetries are non-zero only if CP is violated. Their measurements at future 
colliders allow a determination of the phases, in particular also their signs. 

The asymmetries involve angular distributions of the neutralino, chargino and 
sfermion decay products. The tree-level calculations in the spin-density formal- 
ism include the complete spin correlations between production and decay. Modu- 
lar FORTRAN codes have been programmed for numerical analyses. 
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For neutralino production, e + e — > x° x°j , we can summarize as follows: 

• For the leptonic two-body decay chain of one of the neutralinos x® — > i — > 
X?^2 for ^ = e, /i, r, we have analyzed the asymmetries of two triple products 
7} = (p e - x p^o) • p £l and T H = (p e - x p £a ) • p £l , which are sensitive to and 
</?Mi- In a numerical study for e + e~ — > an d subsequent neutralino decay 
into a right slepton x% — > we have shown that the asymmetry .4/ can be as 
large as 25%. The asymmetry An, which does not require the identification 
of the neutralino momentum, can reach 10%. Asymmetries of the same order 
are obtained for the processes e + e~ — > X1X3 and e + e _ — > x°X3- 

• For the two-body decay of one neutralino into a stau-tau pair, x°i ~^ ^k T , we 
have analyzed the asymmetry of the triple product T T = s T - (p T x p e - ), which 
includes the transverse r polarization s T . The asymmetry is sensitive to the 
phases y? M , Lp Ml and ^a t and can attain values up to 60% for e+e~ — > x?x°- 

• For the neutralino decay into a Z boson, x° — ► -^Xn/ followed by the decay 
Z — > it for £ = e, /i, r, or Z — > gg for q = c,b, we have defined and analyzed 
the asymmetry „4%) of the triple product 7^) = p e - • (p%) x p^), which 
is sensitive to y? M and y? M i • We have also identified the CP sensitive elements 
of the Z spin-density matrix. For X1X2' x\x\> X1X3 and X2X3 production, the 
asymmetry At can go up to 3%. For the hadronic decays of the Z boson, 
larger asymmetries A c (b) — 6.3(4.5) x At are obtained. 

The results for chargino production, e + e~ — > xt xj> can be summarized as follows: 

• For the two-body decay of one of the charginos into a sneutrino xt ~^ 

the asymmetry of the triple product (p e - x p -+ ) • p e is sensitive to y? M and can 
be as large as 30% for e + e" — > xt X2 ■ 

• For the decay xt —> W + x° n the triple product 7} = p e - • (p~+ x p w ) defines 
the asymmetry Ai, which is sensitive to y? M . The asymmetry Ai can reach 
4% for e + e~ — > xt X2 production. Further, we have analyzed the CP sensi- 
tive elements of the W boson spin-density matrix. The phase (p Ml enters in 
the decay xt —> W + x„ due to correlations of the chargino and the W boson 
polarizations, which can be probed via the hadronic decay W + — > cs. More- 
over the triple product T n = p e ■ (p c x p 5 ) defines the cp^ and Lpu x sensitive 
asymmetry An, which can be as large as 7% for X1X1 or X1X2 production. 
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We note in addition that if the neutralinos or charginos decay into scalar particles, 
like sleptons or sneutrinos, the asymmetries probe CP violation in the production 
process only. The asymmetries are then caused by the neutralino or chargino po- 
larizations perpendicular to the production plane, which are non-vanishing only 
for the production of a non-diagonal pair of neutralinos e + e~ — > xl X°j^ 3 or 
charginos e + e — > xf X.2- If the- neutralinos and charginos decay into particles 
with spin, like r or W and Z bosons, we have found that also the CP phases which 
enter in the decay can be probed, in addition to the CP contributions from the 
production. This is due to the spin correlations between production and decay. 

For the two-body decay chain of a sfermion / — > fxj — > fXi Z ~^ f'Xi 1 1 (f'Xi Q 0) 
we have obtained the following results: 

• We have defined the asymmetries At, q of the triple products of the momenta 
of the outgoing leptons % — pf • (p e x p^) or quarks T q — pf • (p q x p^), which 
are sensitive to ip^, and (p Ml - In a numerical study of stau decay t x — > tx\ ff 
we found that the asymmetry At can be of the order of 3% percent for lep- 
tonic final states. The number of produced f's necessary to observe At is at 
least of the order 10 5 , which may be accessible at future collider experiments. 
For the semi-leptonic final state fi — > r Xibb the asymmetry Ab is larger by 
a factor 6.3 and could be measured if the number of produced f's is of the 
order 10 3 . 

• The asymmetries are not sensitive to the phase (pA f , since the decay / — > / x® 
is a two-body decay of a scalar particle. As an observable sensitive to (pA f , 
one would have to measure the transverse polarization of the fermion /, or 
study asymmetries in three-body decays of the sfermion /. 



5.2 Conclusions 

We have shown that in all the processes studied, triple-product asymmetries can 
be defined which are sensitive to the CP phases. Especially promising for mea- 
suring the phases are leptonic decays of charginos and neutralinos, where the 
asymmetries can be as large as 30%. For neutralino decays into a stau-tau pair, 
the asymmetry of the r polarization may even reach 60%. The asymmetries of 
neutralino and chargino decays into a Z or W boson may be as large as 18% and 
7%, respectively. For the neutralino and chargino production processes, we have 
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found that longitudinal polarization of both beams can significantly enhance the 
asymmetries and cross sections. 

By analyzing the statistical significances, we have shown that the asymmetries are 
accessible in future electron-positron linear collider experiments in the 500 GeV - 
800 GeV range with high luminosity and polarized beams. The MSSM phases y? M , 
y?Mi and </?a t could thus be strongly constrained. 



5.3 Outlook 



In further investigations our analysis of CP-odd asymmetries in the production 
and decay of neutralinos and charginos should be extended to electron-positron 
collisions with transverse beam polarizations. For chargino production the asym- 
metries of triple products including the transverse beam polarization are only non- 
vanishing if the decay of both particles is included. 

Moreover, triple products which involve the momenta of the decay products of 
both neutralinos (or charginos) could be sensitive to correlations between their 
polarizations. These spin-spin correlations would yield additional information on 
the CP phases. 



Appendix A 



Basics of the Minimal 
Supersymmetric Standard Model 



The Minimal Supersymmetric Standard Model (MSSM) is characterized by the fol- 
lowing properties: 

• A minimal gauge group: SU(S) x SU{2) x U(l). 

• A minimal content of fields: 3 generations of leptons and quarks, 12 gauge 
bosons, 2 Higgs doublets, and their super partners, see Table A.l. 

• Explicit SUSY breaking parametrized by soft breaking terms. 

• R-parity conservation. 

In the following we give a short account on the relevant Lagrangians and cou- 
plings, and mixing of neutralinos, charginos and staus. 

In terms of superfields, the field content and the parameter content of the MSSM 
is well visible, and the Lagrangian can be written in an elegant and short way. We 
thus give its electroweak part in terms of superfields in Section A.l, where we also 
define the parameters. 

For calculations of interaction amplitudes, however, it is more convenient to ex- 
pand the MSSM Lagrangian in component fields, which is albeit a more complex 
procedure with a longish result. We thus restrict our discussions on the neutralino, 
chargino and stau sector, and give the mass matrices in Section A.3. In Section A.2 
we give the relevant parts of the Lagrangian. 

For detailed reviews of the MSSM and its Lagrangians, see e.g., [58-60]. For a 
detailed study of CP violating sources in the MSSM, see [61]. 
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Table A.l: Particle spectrum of the MSSM 



SM particles 


SUSY-partners 


weak eigenstates 


mass eigenstates 


q u = u,c,t 
q d = d,s,b 


c[lAr squarks 


qi , q 2 squarks 


£ — e, fj,,r 

V = z/ e , V T 


£l,£r sleptons 
i>£ sneutrinos 


l\ , £ 2 sleptons 
vi sneutrinos 


9 


g gluino 


g gluino 


W ± 
(Ht,H 2 ) 


W ± wino 
Hi , H 2 Higgsinos 


xf,2 charginos 


7 
Z° 

Hi (H°) 


7 photino 
Z° zino 
H°, H% Higgsinos 


Xi,...,4 neutralinos 



A.l MSSM Lagrangian in terms of superfields 

We give the electroweak (EW) part of the MSSM Lagrangian in the superfield for- 
malism, which consists of a supersymmetric part and terms which break SUSY 
softly 

Lew = £susy + £ so ft- (A.l) 



A.l.l Superfield content 

The left-handed lepton superfields are arranged in SU (2) doublets and the right- 
handed ones in SU(2) singlets, 

L -(t'S) L - *='*<*•»•»>• < a - 2 > 

where the generation indices of e, \x and r have been suppressed. The super- 
fields are functions on the superspace, with spacetime coordinates and anti- 
commuting Grassmann variables 9 a . The lepton superfield contains both bosonic 



A. 1 MSSM Lagrangian in terms of superfields 
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and fermionic degrees of freedom. It can be expanded 

L = L(x) + iOo^Od^ L(x) - ^OOOWd^ L(x) + 

+V26 L^ix) + -5=000^ L {2 \x) + 00 F L (x), (A3) 

in terms of the spin-0 slepton field L, the spin- 1/2 Weyl field iP^ and a spin-0 aux- 
iliary field F L , which however can be removed by the Euler-Lagrange equations. 
The right-handed field R is expanded similarly to L. 

There are two doublets of chiral Higgs superfields 

where the upper index denotes the SU (2) index. The component field expansions 
of the Higgs fields is similar to that of the lepton field L (A.3). 

The U (1) and 577(2) gauge vector superfields are respectively given by 

V' = Yv'(x,e,e), V = T a V a {x,6,6), (A.5) 

with sum over a = 1, 2, 3 and Y and T a are the U(l) and SU(2) generators. The 
gauge vector superfield contains a bosonic (spin 1) gauge field V®, and a fermionic 
(spin 1/2) gaugino Weyl field A a 

V a (x, 0,0) = - Oa^O V^{x) + iOOO \ a (x) - iOOO X a (x) + -OOOO D a {x). (A.6) 

The auxiliary spin- 1/2 field D a can be removed by the Euler-Lagrange equations. 
The gauge field v' is expanded similarly to V a . 



A.1.2 The supersymmetric Lagrangian 

The supersymmetric Lagrangian contains the kinetic terms for chiral and vector 
superfields, as well as a Higgs part 

£-SUSY — £ Lepton + ^Gauge + £ Higgs- 

(A.7) 
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The lepton and gauge parts are given by 

C Lepton = J d*9 p <?*Wl + # (A.8) 
Ccauge = \ J d 4 9 [W^W* + W' a W' a ] 5 2 (9) + h.c, (A.9) 

where g and g' are the SU (2) and U(l) gauge couplings. The SU (2) and U(l) field 
strengths are defined by 

W a = -^-DDe- 2 ^D a e 2 ^, W' a = --DDD a V', (A.10) 
$g 4 

where D are covariant derivatives. The Higgs part 

Cm 99 s = J d 4 9 [tfje 2 ^ + g'V'H 1 + H\e 2 ^ + g >v'H 2 + W5 2 {9) + W5 2 (9)] (A.11) 
contains the superpotential 

W = W H + W Y . (A.12) 
The Higgs (H) and Yukawa (Y) parts are 

W H = f i &H\Hi W Y = &\JH\UR + hH[Q j D + fcH&U], (A.13) 

with /i the Higgsino mass parameter, f { the Yukawa couplings, with the generation 
index suppressed, and the antisymmetric tensor e 11 = e 22 = 0, e 12 = — e 21 = 1. 
Note that in order to be renormalizable, the superpotential can only be cubic in 
the superfields at its maximum. 

By construction, the Lagrangians given in this section are gauge invariant and in- 
variant under supersymmetry transformations. In addition, they are R-parity con- 
serving, if trilinear terms in W^|e=o are disregarded. The Higgsino mass parameter 
/j, and the Yukawa couplings / can have physical CP-phases. 
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A.1.3 The soft SUSY breaking Lagrangian 



The most general Lagrangian, which breaks SUSY softly, i.e., which does not lead 
to quadratical divergencies, can be classified in scalar mass terms (S) and gaugino 
mass terms (G) [62]: 



C 



soft 



£-S + l~-G- 



(A.14) 



Note that trilinear scalar interaction terms are not R-invariant and will thus be 
neglected. The scalar mass term reads (neglecting squark fields) 



Cs = 



d 4 e 



M\l)L + m\T$R + m\H[H x + m\E\E 2 + 



+mle ij {H\H 3 2 + h.c.)] 5 4 (9, 9), (A.15) 
with the abbreviation M\l)L = m 1 ^ + m 2 L P L l L . The gaugino mass term 

C G = - J d 4 9 \^M x W' a W' a + M 2 W aa W^) + h.c] 5 4 (9, 9) (A.16) 

includes the U(l) and SU(2) gaugino mass parameters Mi and M 2 , respectively. 
By redefining the fields, one parameter, usually M 2 , can be made real, while Mi 
can have a CP violating phase. 



A.2 MSSM Lagrangian in component fields 

The MSSM Lagrangians relevant for chargino production and decay are [58]: 

£zou = 9 —Z^[LzP L + R t P R ]£, (A.17) 

cos 9w 

C ix$x- = -eA^trxt^, e > 0, (A.18) 

= — V^.V[OJPl + OgPdxf, (A.19) 

x-j *» cos 9w J J 

C w-xixl = 9 W ;xlY\0 L kl P L + 0«P R ]xt + h.c, (A.20) 

Ctoxt = -aUaxtPL^l - gV&fPdv* + h.c, t = e, //, (A.21) 

= -glfiV^PL - Y T U l2 P R )ru* + h.c, (A.22) 
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Table A.2: Electric charge e e and weak isospin T 3£ of fermions 





d u e L e R v 




-i +1 -i -i o 

_I +1 _I o 1 

2 ^2 2 U 2 



and the terms relevant also for neutralino production and decay are [58]: 



Zzo^o = l-^-Z^lno'lPL + OZPr]xI (A.23) 

Z, COS u\y 

£ux° = 9&PrX% + 9f^P L X% + h.c, (A.24) 

4,5*2 = gf^ k uP R x k i> L + h.c, (A.25) 

with the couplings 

= T 3 ^ - e^sin 2 ^, = -e^ sin 2 (A.26) 

C# = -l^/i - \v i2 V* 2 + ^ sin 2 (A.27) 

O'if = -U* x U 3l - l -U* 2 U j2 + % sin 2 W , (A.28) 

L kl = -l/V2(cos (3N k4 - sm(3N k3 )V* + (sm6 w N kl + cose w N k2 )V*, (A.29) 

= +l/v^(sin/3JV^ + cosPN* 3 )u i2 + (sm6 w N* kl + cos9 w Nl 2 )u il: (A.30) 

C>™ = -^(iV m3 iV:3 - AUiV: 4 ) cos 2/3 - ^(iV m3 iV* 4 + N m4 N* n3 ) sin2/3, (A.31) 

02 = -O'Z, (A.32) 

^ L fc = -V2\—^—(T 3e - e* sin 2 ^ w )iV fc2 + e £ sin , (A.33) 

.COSC7^ 

fg = -V2e e sin 6 W [ tan %7V* 2 - N* kl ] , (A.34) 

/, L fc = ~V2^—T 3l/ N k2 , (A.35) 



with «, j — 1,2 and k,m,n = 1, . . . , 4. The charge and the third component of the 
weak isospin T 3l of each fermion is given in Tab. (A.2). The r- Yukawa coupling is 
given by Y T — m T / (y/2m w cos (3). 



A.3 Mass matrices 
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For the neutralino decay into staus Xi —> r fc r / we take stau mixing into account 
and write for the Lagrangian [63]: 

C T f Xi = 9hf(al t P R + bl t P L )x° + h.c. , k = 1, 2; % = 1, . . . , 4, (A.36) 
with 

4, = (n f kn )*A T jn , b% = {nl n )*B T jn , 4. = (^L)*oj n , (n = l, a), (A.37) 

^-(5)' 

= (h*)* = -Y T (N* 3 cos (3 + iV* 4 sin /3), (A.39) 

Y T = m T /(V2m w cos(3), (A AO) 

with K f kn given in (A.50) and fff given in (A.33), (A.34). 



A.3 Mass matrices 



A.3.1 Neutralino mass matrix 



The complex symmetric mass matrix of the neutral gauginos and Higgsinos in the 
basis (7, Z, H®, H®) is given by 



Y = 



( M 2 sin 2 6 W + M x cos 2 9 W (M 2 - Mi) sin W cos 9 W \ 

(M 2 - Mi) sin W cos 6> w M 2 cos 2 6> w + Mi sin 2 9 W m z 

t?l2 

V o o 



/i sin 2/5 —fj, cos 2/5 
-/i cos 2/3 — /x sin 2/3 / 



(A.41) 



with /i the Higgsino mass parameter and tan (3 — ^, where u lj2 are the vacuum 
expectation values of the two neutral Higgs fields. The mass matrix Y a/3 (A.41) can 
be diagonalized by a complex, unitary 4x4 matrix [58], 

N* a Y aP N\ k = m t 5 lk , (A.42) 
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with the neutralino masses m^q > 0. Then the weak eigenstates (7, Z, if °, H®) mix 
to the neutralino mass eigenstates (x§, X2, X°, xi)- 

The diagonalization of the neutralino matrix is achieved with the singular value 
decomposition [61]. Let Z be a complex n x n matrix, then: 

• The matrices ZZ^ and Z^Z are selfadjoint and have the same real eigenvalues 

Xi > 0. 

• The eigenvectors e { connected to the eigenvalues A, built an orthonormal sys- 
tem. If some of the A; = 0, the eigenvectors can be completed to built an 
orthonormal system. 

• Let ii be an eigenvector of Z^Z with eigenvalue Aj 7^ 0, that is Z^Z e, = Xi e ir 
then the vectors e- := -^Z also built an orthonormal system with Z^ = 

\fX~iti. 

Thus Z can be decomposed into its singular values: Z — {e[. . . e' 2 )^Z(ei . . . e 2 ) = 
diag(y/Xl). 



A.3.2 Chargino mass matrix 

The complex chargino mass matrix is given by 

X_( M 2 m w \/2sm[3 
cos (3 [i 

It can be diagonalized by two complex unitary 2x2 matrices U mn and V mn [58], 

U* ma X aP V^ 1 = m-+ 5 mn , (A.44) 

with the chargino masses m~+ > 0. The chargino-mass eigenstates xt — (^-) are 

defined by xt = Vnw + + V i2 h + and xj = Uj\w~ + U j2 h~ with w ± and /i ± the two- 
component spinor fields of the W-ino and the charged Higgsinos, respectively. We 
diagonalize the chargino mass matrix with the singular value decomposition, see 
Section A.3.1. 



(A.43) 



A3 Mass matrices 
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A.3.3 Stau mass matrix 



The masses and couplings of the r-sleptons follow from the hermitian 2x2t l — t r 
mass matrix [63]: 



£m — (jhi t r) 



with 



M f LL e-^\M?J\ff L 



e i<Pf\M? LR \ 



M? 

T BB 



(A.45) 



a fl L = M[ + (-- + sin 2 9 W ) cos 2(3 m 2 z + m 2 , 
M? - 



T LL L ' \ 2 

M\ - sin 2 9 W cos 2(3 m| + m 2 , 



M? = (M? )* — m T (A T — //*tan/3), 



arg[A T — // tan (3], 



(A.46) 

(A.47) 
(A.48) 

(A.49) 



where A T is the complex trilinear scalar coupling parameter and M~ L , are the 
other soft SUSY-breaking parameters of the f, system. In order to reduce the num- 
ber of MSSM parameters, we will often use the renormalization group equations 
I + 0.79M 2 and M| = m 2 + 0.23M|. The f mass eigenstates are 



[64], Mf = m 

(ti, f 2 ) = (f L , tr)?^ , with 



/ cos f sin 9 f \ 



n T = 



(A.50) 



\ — sin 9f e l<Pf cos 9 f ) 
with the mixing angle 



cos 9? 



-\M? I 

1 TLB, 1 



mj 2 + (< - m?y 



, sin#f = 



M? 

TLL 



MU 2 + (<-Ml r f 



(A.51) 



and the mass eigenvalues 



rn 



n,2 



( M L + =f 4m~LL - m RR ? + mm? 



2 12 

TLB. I 



(A.52) 
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A3 A First and second generation sf ermion masses 

The off-diagonal terms of the sfermion mass matrices are proportional to the fermion 
mass. For fermions of the first and second generation, whose masses are small 
compared to SUSY masses, their mass matrices are diagonal to a good approx- 
imation. For these sfermions we will assume the approximate solutions to the 
renormalization group equations [64]: 



where T 3 j is the third component of weak isospin, m is the common scalar mass 
parameter at the GUT-scale, and C(f) depends on the sfermion 




mj + m 2 + C(/)M| ± m% cos 2(3(T 3f - e f sin 2 6 W ). 



(A.53) 



C{t L ) « 0.79, C(e R ) w 0.23, 
C(q L ) w 10.8, C(q R ) w 10.1. 



(A.54) 
(A.55) 



For the slepton and sneutrino masses we have 




(A.56) 



(A.57) 



(A.58) 



Appendix B 



Kinematics and phase space 



B.l Spherical trigonometry 

For the parametrization of the phase space one often needs the following relations 
from spherical trigonometry. Consider the following spherical triangle with sides 
a, b, c and angles A, B, C. 




The unit vectors of the triangle sides are given by 

ei = (0,0,1), 



(B.l) 
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e<i = (sin c, 0, cose), (B.2) 
e 3 = (sin b cos A, sin b sin A, cos b). (B.3) 

In the following we give formulas relating the sides and the angles of the triangle 
[65,66]: 

• law of sines 

sin a sin b sin c 



(B.4) 



sin A sin B sin C 
law of cosines for sides (cosine theorem) 

cos a = cos b cos c + sin b sin c cos A. (B.5) 

Similar formulae for the other sides may be obtained by cyclical permuta- 
tions. 

law of cosines for angles 

cos A = — cos B cos C + sin B sin C cos a, (B.6) 
etc, cyclically. 

products of functions from sides and angles 

sin a cos B = cos b sin c — sin b cose cos A, (B.7) 
sin a cos C = cos c sin 6 — sin c cos b cos A. (B.8) 

The products sin b cos C and sin c cos A are obtained from (B.7) by cyclical per- 
mutations. The products sin b cos A and sin c cos 5 are obtained from (B.8) by 
cyclical permutations. 

sin A cos b = cos B sin C + sin £? cos C cos a, (B.9) 
sin A cos c = cos C sin + sin C cos i? cos a. (B.10) 

The products sin B cos c and sin C cos a are obtained from (B.9) by cyclical per- 
mutations. The products sin B cos a and sin C cos b are obtained from (B.10) by 
cyclical permutations. 



B.2 Kinematics of neutralino/chargino production and decay 
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B.2 Kinematics of neutralino/chargino production and 
decay 

B.2.1 Momenta and spin vectors of the production process 

We choose a coordinate frame in the laboratory system (center of mass system) 
such that the momentum of the neutralino x°j or chargino xj, denoted by p Xj , 
points in the ^-direction (in our definitions we follow closely [39,51]). The scat- 
tering angle is 9l(p e - , p Xj ) and the azimuth can be chosen zero. The momenta 
are 

f e _ =E b (l,- sm9, 0, cos 9), p^+ = E b (l, sin0, 0, - cos0), (B.ll) 
p£ = (E Xi , 0, 0, -q), p*. = (E Xj ,0, 0, q), (B.12) 

with the beam energy E b — \f~sj1 and 



s + m 2 — m 2 s + ml — m 2 . A 2 (s, m 2 ., m 2 .) 

Xi Xj_ £j _ Xj Xt_ v ' Xi> Xj ' 



where \(x, y, z) = x 2 +y 2 +z 2 —2(xy+xz+yz). For the description of the polarization 
of the neutralino Xi ° r chargino xt we choose three spin vectors 

4f = (0,-1,0,0), = (0,0,1,0), 4f = JL( qM -E Xi ). (B.14) 

m Xi 

Together with p^./m Xi they form an orthonormal set 

= ^"^=0. (B.15) 



B.2.2 Momenta and spin vectors of leptonic decays 

If the neutralino or chargino decays into a lepton, x° — > t x t or xt —> h&e, in short 
Xi —> hi, it is suitable to parametrize in terms of the angle 6 X = Z(p^, p Xi ): 

p£ = (E £l , - 1 I sin 61 cos 0i, I p ei I sin Q x sin 0i, - 1 pn x | cos 6>i ), (B.16) 
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For the subsequent slepton decay £ — > £ 2 Xi we define 6 2 = l(pe 2 , p Xi ) and write 

p£ = (£^ 2 ,-|p/ 2 |siii02COS02, |p/ 2 |sin0 2 sin02,-|p< 2 |cos02), (B.18) 

Ee 2 = Pe 2 = ^= - 1 * (B.19) 

2{Ef- \p Xl - p ei \cos9 D2 ) 

with 6 2 = l(pe 2 , p Xi ), the decay angles 9 D2 L{p b p ta ), 9 Dl l(p Xi , pi) and 

cosfl D2 = cos# Dl cos^2-sin^ Dl sin 6 2 cos(0 2 -0i), cos6 Dl = Pxi [ Px ' ~ (B.20) 

IPxillPxi ~P<il 

If the neutralino decays into a stau, x° — ► ^i", m = 1,2, the r spin vectors are 
chosen by 

1 / S 2 X S 3 \ 2 / Pr X p e - \ 3 I ( E T \ 

S r= °> I 3 1 . S r = 0,, 1 , 4 = — ( |Pt|,— Pr I ■ ( B " 21 ) 



S 2 X S3 J V p T x p e -| / m T V |p 



T I 



B.2.3 Phase space for leptonic decays 

For neutralino /chargino production e+e — > XiXj an d subsequent leptonic decay 
Xi —> £\£ or xt ~> h^e, i n short Xi —> hi, the Lorentz invariant phase-space ele- 
ment can be decomposed into two-body phase-space elements [65]: 

dLips(s; p x . , p tl , p^) = — dLips(s; p Xi , p x . ) ds Xi dLips(s Xi ; , p^) . (B.22) 

For i = I, we can include the subsequent selectron decay £ — > £ 2 X? and have for 
the complete process e + e~ — > X° - *• hh £ —> £2X1'- 

dUps(s;p Xj ,p h ,p i2 ,p x o) = — dUps(s;p Xj ,p ei ,pi) ds g dUps(sfpe 2 ,p x o,). (B.23) 

The several parts of the phase space elements are 

q 

dUps(s;p Xi ,p Xj ) = sin0 dO, (B.24) 

1 \Pi I 2 

dhips(s Xi ;p^p ei ) = — jdfii, (B.25) 
s 2{2ny m x . — m~ 

1 |p^ 1 2 

dLips(s^ 2 ,p 0) = . . 2 2 rifl 2 , (B.26) 
Al 2(2n) 2 mj — m Xl 



B.2 Kinematics of neutralino/chargino production and decay 
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with s Xi = p x _, = p~ and dVti = sin#; ddi dfa. We use the narrow width ap- 
proximation for the propagators / |A(xi)| 2 ds Xi = m * v ^ , / |A(f)| 2 c^~ = The 

approximation is justified for (T x Jm Xi ) 2 1, and (I^/m^) 2 <C 1, which holds in 
our case with T Xi < C(lGeV) and T i < C(lGeV). 



B.2.4 Energy distributions of the decay leptons 

For neutralino production e + e~ — > and the subsequent leptonic decay x° —> 
£i£, and £ — > £2X1/ the two decay leptons £1 and £ 2 can be distinguished by their dif- 
ferent energy distributions. The energy distribution of lepton i x in the laboratory 
system has the form of a box with the endpoints 



2 2 
m — to- 

Eil,min,max <-> / n i \> (B.27) 

2(E x o ±g) 



with q the neutralino momentum. The energy distribution of the second lepton £ 2 
is obtained by integrating over the energy E~ k of the decaying slepton 



, > I Zn^r 5 A <E h < a 

1 drr mi TO 4 ~~ ~~ 

' " 1 xl - J ln\ ; a <E i2 < b (B.28) 



a d£ £2 g[m 2 - to?] [to? - m\] 



with 



Zn-r?- ; 6 < Ef < B 



2 2 

TO- — TO ~,o 



A ^ = (^ >fnax =f V^- m f) (B - 29) 



2to| 
to? - m\ 



E x o{m\ + to?) ± (m 2 - to?) JE\ - m\ 



^i^ax.min ~ 2TO 2 ' \D.OL) 

X? 
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B.2.5 Momenta and spin vectors of bosonic decays 

For the bosonic two-body decays of neutralino x° — > Z°Xn or chargino xt ~> 
W + Xn> i R short Xi —> Bx Q n , we define the decay angle between neutralino (or 
chargino) and the boson B as 6 1 l(p Xi , ps). The angle is constrained by sin fl™** = 
g'/g for g > g', where g' = A 2 (m 2 . , m^, m 2 )/2ms is the neutralino (chargino) mo- 
mentum if the boson 5 is produced at rest. In this case there are two solutions 



± (to 2 . + m 2 B - m x o)qcos9i ± E Xi yjX(m x .,m%,m x o) - ^q 2 m 2 B (1 - cos 2 #1) 
' Ps ' = 2g 2 (l - cos 2 6 1 ) + 2m x . ' 

(B.32) 

For q' > q, the angle Q\ is not constrained and only the physical solution |p#| is left. 
We parametrize the momenta of the decay B — > // in the laboratory system as 

p^ ,lx = (E B , — \p B \ sin cos 0i, \p B \ sin^iSin^!, — |p^| cos#i), (B.33) 
= (£■/,- 1 p/l sin 6*2 cos 02, |p/| sin6> 2 sin02, -|p/| cos6> 2 ), (B.34) 

with 9 2 /.(p Xi ,Pf) and the decay angle 6 D2 l(p B , pj) given by 

cos 9 D . 2 = cos #1 cos 6 2 + sin #i sin 9 2 cos(0 2 — 0i) ■ (B.36) 
The three spin vectors of the boson B = Z°, W + are in the laboratory system 



V ItflXtsI/ V |Pe-xp B |y to s V |p B | 



(B.37) 



Together with p B /rriB they form an orthonormal set. The polarization four- vectors 
e Afc for helicities A fc = —1, 0, +1 of the boson are defined by 

£ " = 73(^-^1); e = f|; e + = -^ + t^). (B.38) 



B.3 Kinematics of sfermion decays 
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B.2.6 Phase space for bosonic decays 



For neutralino/chargino production e + e~ — > XiXj an d subsequent decay of the 
neutralino x° — ► or chargino x+ — > W + x®, in short x« — > #X°' tne Lorentz 

invariant phase-space element can be decomposed into two-body phase-space el- 
ements [65]: 

dUps(s;p Xj ,p x o n ,p B ) = ^- dLips(s;p Xi ,p x .) ds Xi ^ dLips(s Xi ;p x o (B.39) 



If we include the subsequent decay B — > // we have 

rfLips(s;p Xj ,p x o,p / ,p/) = rfLips(s;p Xj ,p x o,p B ) ds B rfLips(s B ;p / ,p / -). (B.40) 

The several parts of the phase space elements are given by 

q 

dUps(s;p Xt ,p Xj ) = ^— y= sm9 d9, (B.41) 

dU P s(s x -p x o,p±) = ^ 2W ± q J*f_ E ^ m (B.42) 

1 |p,-| 2 

dUps(s B ;p f ,Pf) = dft 2 , (B.43) 

with s Xi = p Xi , sb = Pb an d dfli = sin0j dfa. We use the narrow width ap- 
proximation for the propagators / |A(x;)| 2 ds Xi = ^ JI r ~, J \A(B)\ 2 ds B = 

Xi Xi ts t> 

The approximation is justified for (Y Xi /m Xi ) 2 1, which holds in our case with 

r Xi <G(lGeV). 



B.3 Kinematics of sfermion decays 
B.3.1 Momenta and spin vectors 

We consider the slepton decay chain t — > Ix®, x® ~^ Xi^, Z — > //. The substitu- 
tions which must be made for similar decay chains of a squark are obvious. The 
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momenta in the slepton rest frame are 

Pz = (£ z ,0,0,|p||), (B.44) 

Pto = \Pyo\(E Y g/|p Y q|,sin0i,O,cos0i), (B.45) 

Pf = | P/K-E 1 // 1 P/l, sin ^ 2 cos 02, sin ^ 2 sin 02, cos ^ 2 ), (B.46) 
2 2 

|P *°' = -^T' |P/I = 2(E z -|p z |cose 2 )- (B - 47) 

There are two solutions for |pf |, see (B.32), if the decay angle Q\ = Z(pz,p Xi ) is 
constrained by 

m,- X^(m1,m 2 7: m 2 , ) 

sin ^[ nax = — 7 ¥ o ^ < 1 • (B.48) 

™z {mj-m x .) 



The y° spin vectors in the £ rest frame are 



Together withp^p/^ x p they form an orthonormal set. 



(B.49) 



B.3.2 Phase space for sfermion decays 

The Lorentz invariant phase-space element for the decay chain £ — > £x° jr x° — > 
Z —> ff can be written in the rest frame of £ as 

dLips(mf; p*, p x? , p 7 , p/) = 

j^y dUps(mf,p e ,p x o) ds x o^2dUps(s x o;p x o,p^) ds z dUps(s z ;Pf,p f ), (B.50) 



B.3 Kinematics of sfermion decays 
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1 / . m l° 



8(2vr) 2 l m 





Pt! 


2 






COS^i 




Pill 



dLips(m|;p € ,p x o) = ( 1 - -3/ I dQ > ( B - 51 ) 
dLips(s x o;p x o,p|) = |p x0 '-^' ^ i-±n (B " 52) 

dLips^zip^p/) = ± j-J -^ d ^2, (B.53) 

J 8(27r) 2 (Ez — \Pz\ cos0 2 ) 

with s o = p 2 , s z = p% and dVti = sin#; dOi d<pi. We use the narrow width ap- 
proximation for the propagators / |A(x°)| 2 ds x o = ^-jf^, J \A(Z)\ 2 ds z = 

The approximation is justified for (T x o/m x o) 2 <C 1, which holds in our case with 

r x o<G(lGeV). 



Appendix C 

Spin-density matrices for neutralino 

production and decay 



We give the analytic formulae for the squared amplitudes for neutralino produc- 
tion, e + e~ — > x° x°/ with longitudinally polarized beams, and for different sub- 
sequent two-body decay chains of one neutralino. We use the spin density matrix 
formalism as in [27, 39, 66]. The amplitude squared can then be written 

|T| 2 = |A(x°)| 2 £ Pp(x^ K Pn{x°) K x v (CI) 

with pp(x°) the spin density production matrix of neutralino x°, the propagator 

A(x°) = i/[s x q — m 2 x o + im x oT x o] and the neutralino decay matrix Pd{xT)- 



C.l Neutralino production 

For the production of neutralinos 

e + + e - x°(p x o, Aj) + x?(p x o, A,), (C.2) 

with momentum p and helicity A, the unnormalized spin-density matrix of neu- 
tralino Xi is defined as 

Pp(x°) X ^ = J2Tp iXj Tp iXj *. (C.3) 



C.l Neutralino production 
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The helicity amplitudes are [39, 66]: 



Tp lAj (s, Z) 



T* iXj (t,e L ) 
T* iXj (t,e R ) 



Tp' X \u,e L ) 



Tp iX] (u, e R ) 



-^L—A s (Z)v( Pe+ )Y(L e PL + R e P R )u{p e -) 

COS^ Vw 

xu(p x o, X^O'^Pl + 0"^P R )v{p xV \), 

-9 2 feifej ¥At (e L )v(p e +)P R v(p x 0,X i )u(p x 0, \j)P L u(p e -), 

-g 2 f?J**k\e R )v(p e+ )P L v(p xV X l )u(p x o, \j)P R u(p e -), 
g 2 f^f^ u {e L )v{p e+ )P R v{p xV \j)u{p^, \i)P L u(p e -), 
g 2 f^f^ u {e R )v{p e+ )P L v{p x ^ A^, Xi)P R u(p e -). 



(C.4) 
(C.5) 
(C.6) 
(C.7) 
(C.8) 



The propagators are 



A S (Z) = 



s — m| + imz^z 



, AW) 



, AW) 



-, (C.9) 



with s = (p e - +p e +) 2 , t = (p e - —p x o) 2 and u = (p e - —p x o) 2 . The Feynman diagrams 
are shown in Fig. C.l. 







e L,R 




Figure C.l: Feynman diagrams for neutralino production 



For the polarization of the decaying neutralino x° with momentum p x o we have 
introduced three space like spin vectors s a (B.14). Then the neutralino production 
matrix (C.3) can be expanded in terms of the Pauli matrices, see Appendix F.1: 



0\A;A' 



2(5 AiA ,P + ^,S») ) 



(CIO) 



where we sum over a and the factor 2 is due to the summation of the helicities of 
the second neutralino x°, whose decay will not be considered. With our choice of 
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the spin vectors, £p/P is the longitudinal polarization of neutralino x°, £p/P is 
the transverse polarization in the production plane and Sp/P is the polarization 
perpendicular to the production plane. Only if there are non-vanishing CP phases 
Lp Ml and/ or ip^ in the neutralino sector, and only if two different neutralinos are 
produced, e + e~ — > Xi Xj/ i 7^ h t ne polarization Ep/P perpendicular to the pro- 
duction plane is non-zero. Thus it is a probe for CP violation in the production of 
an unequal pair of neutralinos. Note that Sp also gets contributions from the finite 
Z width, which however do not signal CP violation. 

We give the analytical formulae for P and Sp, Sp, Sp in the laboratory system [66] 
in the following sections. Lorentz invariant expressions for these functions can be 
found in [39, 66]. 



C.l.l Neutralino polarization independent quantities 

The coefficient P is independent of the neutralino polarization. It can be composed 
into contributions from the different production channels 



P = P(zTzT) + P{Z~e R ) + P{e R e R ) + P{Z~e L ) + P(e L e L ), 



(C.ll) 



with 

P(zTzT) 



9 



COS* V W 



A'(Z)| 2 W* + L&03 



x{\Ol R \ 2 (E x oE x o + q 2 cos 2 9) - [{ReO'^f - (ImO'fflm^m^}, (C.12) 



P(Ze R ) = -JL^R eCR E 2 b Re{A s (Z) 



COS^ U w 



X 



R rR*/~\"R\ 

ei J ej ij 



+(A*(g fl )/*/*0" i R + A w (e R )f£fS> 0':«* )(E x oE x o + q 2 cos 2 6) 



R fR*r^"R*^ 



eUej ij 



-(A**(e R )/ e f /*oy - A"(e fl )/5/*O^)2^008fl]}, 

= ^c«^{|/5| 2 |/5l 2 >< 

;(|A*(ep)| 2 + \A»(e R )\ 2 )(E x oE x o + q 2 cos 2 9) 

-(\A\e R )\ 2 -\A u (e R )\ 2 )2E bq cos9 
-Pe{(/ e f) 2 (/5) 2 A"(e R )A t *(ep)}2m x om x o}. 



(C.13) 



(C.14) 



C.l Neutralino production 
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To obtain the quantities P(Ze L ), P{e L e L ) one has to exchange in (C.13) and (C.14) 



A\e R ) - A\e L ), A u (e R ) - A u (e L ), c R - c L 

p , T f\"R __. f)"L f R f L f R cL 

■"-e -^e; ^jj ^ij ; ,7 e j ,7 e j; j e j Jej- 



(C.15) 



The longitudinal beam polarizations are included in the weighting factors 

C L =(1-P c -)(1 + P e +), Cfl=(l+P e -)(1-P e +)- ( C16 ) 

Generally the contributions from the exchange of e R (e£) is enhanced and that of 
e L (e R ) is suppressed for P e - > 0, P e + < (P e - < 0, P e + > 0). 



C.1.2 Neutralino polarization 

The coefficients E P , which describe the polarization of the neutralino decom- 
pose into 



E P = Ep(ZZ) + Ep(Zgp) + Ep(gpgp) + E P (Zg L ) + E P (g L g L ). 



(C.17) 



• The contributions to the transverse polarization in the production plane are 



Y} P (ZZ) = 2 



9 



cos 4 9 W 



A s (Z)\ 2 E 2 sm9(R 2 e c R -L 2 e c L ) 



x 



Ep(Zgp) = 



Ol R \ 2 m x oE x , - [{ReO'^f - (ImO'^m^ 
~9 4 



(C.18) 



-R e c R E b sin 9 



x 



- Re{A s (Z)[ffjf0:f*A u *(i R ) + f** /«O^A**(ep)]m x o£; x o} 

+Re{A s (Z)[f e R J e fO'; j R A u *(e R ) + f^A^^m^}], (C.19) 
5 



SM^) = ycp^sin^{|^| 2 |/e1l 



(|A'(gp)| 2 + |A u (gp)| 2 )m x oE x o - (|A*(gp)| 2 - |A u (ep)| 2 )m x o ? cos0 



-2 J Re{(/ e t) 2 (/5) 2 A«(gp)A'*(gp)}m x oE x o}. 



(C.20) 
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To obtain the expressions for T} p (Ze L ) and T} p (e L e L ) one has to apply the ex- 
changes (C.15) in (C.19) and (C.20) and to change the overall sign of the right hand 
side of (C.19) and (C.20). 

• The contributions to the transverse xl polarization perpendicular to the pro- 
duction plane are 



E 2 P (ZZ) 
Z 2 P (Ze R ) 



cos 4 6w 
A 



\A s {Z)\ 2 {R 2 e c R - L 2 e c L )m x oqE 2 S m9Re(0" R )Im(0'' R ), (C.21) 



, -R e c R m x qE 2 q sin 9 
cos 2 9w 3 

x Im{A s (Z)[f*f e f 0'**r(e R ) - / e f / e ^*A'*(e R )]}, 



X 2 P (e R e R ) = -|c R m x o^gsin^/m{(/ e t) 2 (/5) 2 A«(g i? )A t *(g R )}. 



(C.22) 
(C.23) 



To obtain the expressions for Tip(Ze L ) and Tip(e L e L ) one has to apply the ex- 
changes (C.15) in (C.22) and (C.23). 

• The contributions to the longitudinal Xa polarization are 



Y? P {ZZ) 



T? P {Z~e R ) 



2g 4 



cos 4 9 W 



\A s (Z)\ 2 (L 2 e c L -R 2 e c R )E 2 b cos9 



x 



OlfnE x oE x o + - [{ReO" R f - (ImO;fr]rn x om x « 

4 



V'R*\2-i 



(C.24) 



-g 



-R e c R E b 



cos 2 9 W 

x [ite{A-(Z)[/*/*O^A"(e fl ) - ftfp:*^(e R )}E x »q) 
+Re{A%Z)if e R jfO:**A-(e R ) + f R * / e ^A<*(e«)](£ x oi? x o + q 2 ) costf} 
-Re{A'(Z)\f^O^(e R ) + f R 0'^ A^{e R )\m x orn x o costf} 
+Re{A s (Z)[f e R J e fO'; j R *A u *(e R ) - f R * f R 0': R A* (e R )]E x oq cos 2 9}] , (C.25) 



Z P (e R e R ) = 



-^c^[i/5i 2 i/5i 2 

xf[|A«(e fl )| 2 - |A*(g i? )| 2 ]Kog+ [|A«(e~ R )| 2 - \A\e R )\ 2 \qE x » cos 2 



+ [|A'(^)| 2 + |A u (e*)| 2 ](£ x o£ x o + g 2 ) costf} 
-2i?e{(/ e t) 2 (a 2 A«(g i? )A**(e^)}m x om x ocos^ 



(C.26) 
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To obtain the expressions for Y? p {Ze L ) and Y? P {e L e L ) one has to apply the ex- 
changes (C.15) in (C.25) and (C.26) and to change the overall sign of the right hand 
side of (C.25) and (C.26). 



C.2 Neutralino decay into sleptons 

For neutralino two-body decay into sleptons 

Xi(p x o,Xi) - l + t = e,f,,T, (C.27) 
the neutralino decay matrix (2.8) is given by 

PoAx^Kx^h^ + at,^, (C.28) 



where we sum over a. For the decay into right sleptons x° — »■ if, t — e, \i, the 
expansion coefficients are 

Di = y|/^ 2 Ko-mf), (C.29) 
= ±g 2 \f*\ 2 m x o(s a x o- Pei ). (C.30) 

For the decay into the left sleptons xl —> ^f, i = e,fi, the coefficients are 

Di = £|/£lVJo-mf), (C.31) 



2 



^ = T9Vk\ 2 m x o(s xr p ei ). (C.32) 

For the decay into the stau xf ^ ^ T± > k = 1,2, one obtains 

Di = y (14/ + l&LHKo - m^J, (C.33) 

S a Dl = ^ 2 (|4,| 2 - \bU 2 )m x o(s x o- Pei ). (C.34) 

The factor for the subsequent slepton decays — > ^2X1/ i = e,/i, is given by 

D 2 = <f|/£^(ro|-roJo), (C.35) 
and that for stau decay f k — > rx5 by 

^2 = ^ 2 (l4il 2 + l&Ll 2 )K-^ 2 o). (c.36) 
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C.3 Neutralino decay into staus 

For neutralino two-body decay into staus 

X°(p x o, Aj) - f ±(p f J + r T (p T , A fc ); m = 1, 2, (C.37) 
the decay matrix is 

Mx°)$f = S KXl D^ + £^(E^)^. (C.38) 

a 

With the spin basis vectors for the r T , given in (B.21), we can expand 

D x ^ = 5 Xk yD + a b XkK D\ (C.39) 

The expansion coefficients are given by 

D = g 2 Mbi;al^m T m x o + ^(|&^| 2 + |<4| 2 )(p T • p x o), (C.41) 

D b = ± g -^mMj 2 - \aU 2 )(p x o ■ s b T ), (C.42) 

= ±^m x o(|aL| 2 - |&L| 2 )(Pr • s a x o), (C.43) 
^ = /Re(6L*aL)(K-^o)(p x o-4) 

-^ 2 (^o • s b T )[^(\b T mi \ 2 + \a T mi \ 2 )m T m x o + Re(b T m *a T mi ){p T ■ p x o)} 

T g 2 Im(bl*a^)e, upa ptp x oS a /s b T '°, (e 0123 = 1). (C.44) 



C.4 Neutralino decay into the Z boson 

For the neutralino two-body decay into the Z boson 



C.4 Neutralino decay into the Z boson 



113 



the decay matrix is given by 

with the helicity amplitude 

= u(p x o, W-^-P^Pl + OjP R ]u(p x o, X^*. (C.47) 
cob 

For the subsequent decay of the Z boson 

Z(pz, A fc ) - /(p/, A/) + A/); / = £, g, (C.48) 
the decay matrix is 

P D2 (Z) KXk = £ <XT X D % (C.49) 
A /' A / 

and the helicity amplitude 

^ = «(P/> X f )-f-^-[L f P L + ^Mp/, ( C - 5 °) 



The polarization vectors of the Z boson e^ k ,Xk = 0, ±1, are given in (B.38). With 
the set of neutralino spin vectors s%, given in (B.14), we obtain for the neutralino 
decay matrix 

pnAtt$ = o^r + vi^iryr^' ( c - 51 ) 

and for the Z decay matrix 

Pd 2 (Z) x , Xk = D?ej*ei"; (C.52) 
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with 



2g 2 



cos 2 6w 



{[2 p^o - (f xfz + v^fz) ~ IKo + m 2 o - m 2 z ) 9 n\0" n 



L|2 



-g^m x om x0n [(ReOjf - (ImOjf]}, 

{-m x ^sl la (p xlp - Pz , p )\0: 



(C.53) 



2^ 2 



cos 2 # 



L|2 
ni I 



+2m x o - s^p^)(ImO n f)(ReO. 



—m x o e 



navf3 



ni ) 



s a x>Px oA(ReO':f) 2 - (ImO':t) 2 }}; (e i 23 = 1), 



(C.54) 



and 



2^ 2 



{(-2 pjp£- + p^- + pjp£ - |m|<T)(Lj + 4) 



cos 2 W 



(C.55) 



Due to the Majorana character of the neutralinos, is symmetric and Ep^" is 
antisymmetric under interchange of fi and v. In (C.51) and (C.52) we use the ex- 
pansion (E12) for the Z polarization vectors 



2m z 



(C.56) 



summed over c, d. The decay matrices can be expanded in terms of the spin matri- 
ces J c and J cd , given in Appendix R2. The first term of the decay matrix p Dl (C.51), 
which is independent of the neutralino polarization, then gives 



cd\X k X' 



(C.57) 



with 

D 1 



cos 2 9 



w 



An 



m| + 



4 (p x o • Pz) ; 



m 



10 



'L|2 



+2m x om x o [(.ReO^) 2 - (/mOl L ) 2 ]}, 



{ [2(4 " P x? )(4 • p x o) + ±(m 2 „ + m 2 ( , - m|)5 cd 



COS 2 #w 

+5 cd m x om x o [(ReOjy - (JmOlf) 2 ]}, 



|0" L | 2 



(C.58) 



(C.59) 
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and c Di = due to the Majorana character of the neutralinos. As a consequence of 
the completeness relation (F.14), the diagonal coefficients are linearly dependent 



11 



D 1 + 22 D 1 + 6S D 



33 



-ID,. 



(C.60) 



For large three momentum p x o, the Z boson will mainly be emitted into the for- 
ward direction with respect to p x o, i.e. p x o pz p z , with p = p/|p|, so that (t^ 2 -p x o) pz 
in (C.59). Therefore, for high energies U D ± pz 22 D lr and the contributions of the 
non-diagonal coefficients cd Di (c ^ d) will be small. 

For the second term of pr> 1 (C.51), which depends on the polarization of the decay- 
ing neutralino, we obtain 



with 



(C.61) 



2g 2 



m z cos 2 9 W 



{[\OZ \ 2 m x o + i(ReO:t) 2 - (ImOj) 2 ^ 



s a x o ■ pzWz ■ P x o) - (s a x o ■ t%){p z ■ p x o)} + \Oj \ 2 m x om 2 z (s a x o ■ t%) 
-2(/mO:f)(i?eO:f)m x o6, 1 , p(7S ^„/ z t c /} ) (C62) 



and YP D 



cd 



Y> a Di = due to the Majorana character of the neutralinos. Insert- 
ing (C.57) and (C.61) into (C.51), we obtain the expansion of the neutralino decay 
matrix 



(C.63) 



into the scalar (first term), vector (second term) and tensor part (third term). 
A similar expansion for the Z decay matrix (C.52) results in 

pD 2 (Z) K x k = D 2 5 X '^ + C D 2 (J c f* x « + cd D 2 (J *)^*, (C.64) 
with 

2g 2 



D 2 
cD 2 



3 cos 2 9 W 
2g 2 



(Rj + L 2 f )m 2 z , 



cos 2 9w 

g 2 



(Rj-L 2 f )m z (t z -p f ), 



cd n y ( D 2 



cos 2 9w 



(Rj + Lj) 2(t z -p f )(t d -Pf) 



\m\b cd 



(C.65) 
(C.66) 
(C.67) 
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As a consequence of the completeness relation (F.14), the diagonal coefficients are 
linearly dependent 

n D 2 + 22 D 2 + 33 D 2 = -\D 2 . (C.68) 

For large three-momentum pz, the fermion / will mainly be emitted into the for- 
ward direction with respect to pz, i.e. pz ~ p/, so that (t^ -pj) ~ 0in(C.67). There- 
fore, for high energies U D 2 m 22 D 2 , and the contributions for the non-diagonal 
coefficients cd D 2 (c ^ d) will be small. 



Appendix D 

Spin-density matrices for chargino 

production and decay 



We give the analytic formulae for the squared amplitudes for chargino production 
e + e~ — > xt Xjf with longitudinally polarized beams and for different subsequent 
two-body decay chains of one chargino. We use the spin density matrix formalism 
as in [27, 51, 66]. The amplitude squared can then be written 

|T| 2 = \A(xt)\ 2 E Pp(*t) XtK PnWhx, (D-D 

with pp(xt) the spin density production matrix of chargino xt' the propagator 

A(xt) — V[ s v+ — r m 2 ++ im +T +] and the chargino decay matrix Pd{xD- 



D.l Chargino production 

For the production of charginos 

e + + e~ -> (p x + , Ai) + xj (p x j Aj) » ( D - 2 ) 

with momentum p and helicity A, the unnormalized spin-density matrix of chargino 
xt is defined as 

Pp(xt) XlK =Y, T p iXj T^*. (D.3) 



118 



Spin-density matrices for chargino production and decay 



The helicity amplitudes are [51, 66]: 



T^( 7 ) = -e 2 A( 7 )(J -i;(p e+ )y t «(p e -)iZ(p Y +,A i )7^(p Y 7,A J -), 



Tp lX] (Z) = 



-A(Z)v( Pe+ )r(L e P L + R e P R )u{p e -) 



cos 2 9w 

xu(p x +, X t )^(0' t jP L + 0\fP R )v(p x -A 3 ), 
T X p ^(u) = VF i i^A(i/)i;(p e+ )Pfl7;(p x +,A0«(p x7 ,A J -)PL«(Pe-), 



with the propagators 



A(7) = 4» A ( Z ) = - 2' • f 

P 2 Pz~ m z + %m z^z 



Pi -mi 



The Feynman diagrams are shown in Fig. D.l. 



(D.4) 

(D.5) 
(D.6) 



(D.7) 



Xi 



Xi 



X, 



V 



1/ I 



A 



x, 



Figure D.l: Feynman diagrams for chargino production 

For the polarization of the decaying chargino xt with momentum p x + we have 
introduced three space like spin vectors s a + (B.14). Then the chargino production 
matrix (D.l) can be expanded in terms of the Pauli matrices, see Appendix El: 



p P (xt) X ^ = 2(5 AiA ,P + < A ,£ P ), 



(D.8) 



where we sum over a. The factor 2 in (D.8) is due to the summation of the helicities 
of the second chargino xj , whose decay will not be considered. With our choice 
of the spin vectors, Xp/P is the longitudinal polarization of chargino xt> ^p/P 
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is the transverse polarization in the production plane and Sp/P is the polariza- 
tion perpendicular to the production plane. Only if there is a non-vanishing CP 
phase y? M in the chargino sector, and only if two different charginos are produced, 
e + e~ — > xf X.2' t ne polarization Ep/P perpendicular to the production plane is 
non-zero. Thus it is a probe for CP violation in the production of an unequal pair 
of charginos. Note that Ep also gets contributions from the finite Z width, which 
however do not signal CP violation. 

We give the analytical formulae for P and Ep, Ep, Ep in the laboratory system in 
the following sections. Lorentz invariant expressions for these functions can be 
found in [51,66]. 



D.l.l Chargino polarization independent quantities 



The coefficient P is independent of the chargino polarization. It can be composed 
into contributions from the different production channels 



P = p( 77 ) + P( 7 Z) + P(7*>) + P(ZZ) + P{Zu) + P{Vv) 
which read 



(D.9) 



P( 77 ) 
P{lZ) 

P(ZZ) 



P(Z9) 



P{vv) 



<%2e 4 | A( 7 )| 2 (c L + c R )E 2 b (E x +E x -+m xt m xT +q 2 cos 2 9), 



(D.10) 



c e 2 g 2 
<%2 — 2~o 

COS u\y 



P 2 Pe{A( 7 )A(zT)* [(L e c L - R e c R )(0' t f* - 0%*)2E b q cos 9 



+ (L e c L + R e c R )(0' t f* + 0' t f*)(E x+ E x - + m x+ m x - + q 2 cos 2 0)1 }, (D.ll) 

1 J 1 J -Is 

^P 2 c LJ Re{\/; i \/ jl A( 7 )A(z>)*} 



dije 



x(E x +E x - +m Y +m x - - 2E b qcos9 + q 2 cos 2 9), 

4 



9 



(D.12) 



cos* v w 

+(Llc L + Rlc R )(\0'£ 



\A(Z)\ 2 E 2 (L 2 e c L - R 2 e c R )(\0' t f | 2 - \0^\ 2 )2E h qcos9 



l 2 + 


\o' R \ 




1 *3 1 



){E y+ E-+q 2 cos 2 9) 



+ (L 2 e c L + Rlc R )2Re{0^0'«*}m xt m xZ 



-L e c L E 2 Re{v* l V :il A(Z)A( 



cos 2 9 W ~ 

x[0'^{E x +E x - - 2E b qcos9 + q 2 cos 2 9) + O'.fm^m^}} , 



^c L \V^ 2 \VM^)\ 2 El{E xt E x - 



2E b qcos9 + q 2 cos 2 9). 



(D.13) 

(D.14) 
(D.15) 
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The longitudinal beam polarizations are included in the weighting factors 

c L = (1 - P e -)(1 + P e+ ), c R =(l + P e -)(1 - P e+ ). (D.16) 
Sneutrino exchange is enhanced for P e - < and P e + > 0. 

D.1.2 Chargino polarization 

The coefficients £ P/ which describe the polarization of the chargino xt > decompose 
into 

S P = S P ( 77 ) + S P ( 7 Z) + S P ( 7 z>) + Sp(ZZ) + S P (Zz>) + S P (z>z>). (D.17) 
• The contributions to the transverse polarization in the production plane are 



Sp( 77 ) 

Sp( 7 zT) 



Sp( 7 z>) 
T} P (ZZ) 



Y} p {Zv) 



£p(z>£) 



%2e 4 |A( 7 )| 2 ( Ci? - c L )E 2 b sin 0(m x +E x - + m x -E x +), 
^2^^^ 2 sm^e{A( 7 )A(zT)* 



(D.18) 



•'ij 
x 



cos 2 9 W 

- (L e c L + R e c R )(0' t f* - 0'^*)m x +qcose 

E v - +m^-E xt )]}, 



+(i? eCi? - L e c L )(0;f + 0'£*)(m 4 E Xj 
P 2 sin^e{\/; i \/ jl A( 7 )A(z>)*} 



r 2 2 



X 



(D.19) 
(D.20) 



2 _ 






\^IJ 1 



COS 4 ^VK 

+(it! 2 C p_ - L 2 c L )2i?e{0; L 0;f }m x -P Xi+ 
+(i? 2 c* - L 2 c L )(|0;f | 2 + \0'i\ 2 )m xt E x - 

? n -L e c L E 2 b sm6Re{v*V :il A(Z)A(vy 

COS z [7 14/ k 

-^c L |^ri^ 1 | 2 |A(z>)| 2 P 2 sin^m x +(P x - -gcos^). 



30 



(D.21) 

(D.22) 
(D.23) 
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The contributions to the transverse xt polarization perpendicular to the pro- 
duction plane are 



£ 2 P ( 77 ) = Z 2 P (vv) = 0, 



6 ^ 2 ^h {ReCR ~ L eCL)/m{A( 7 )A(zT)*(0;f - O'f )} 



xE 2 m -q sin 9, 



£p(7*>) = 5 ij e 2 g 2 c L Im{V*V jl A( 1 )A(v)*}Etm x -qsin9, 
Ep(zTzT) 



J cos 4 9 W 



\A(Z)\ 2 (R 2 e c R - L 2 c L )/m{0; L 0;f }E 2 m x -gsin^ 



cos 2 9 W 



L e c L Im{v* 1 V ll O' i fA(Z)A(i>y}E 2 b m xT qsm9. 



The contributions to the longitudinal xt polarization are 



(D.24) 

(D.25) 
(D.26) 

(D.27) 
(D.28) 



£2.(77) 

S 3 P ( 7 ZT) 



T?p(ZZ) 



Y? P {ZV) 



Y? P {vv) 



<%2e 4 |A( 7 )| 2 (c L - c p )£ 2 cos% 2 + ^ x +^ x - + m y +m y7 ), 



%2 — — 



S 2 i?e{A( 7 )A(zT)^ 



x 



cos^ v w 

(L e c L - R e c R ){0'£* + 0'^){q 2 + E xt E x - + m xt m x -) cos9 



+ (L e c L + R e c R )(0' t f* - 0' t f*)q(E x - + E xt cos 2 9)} }, 
^pM = -%e 2 /c L E 2 J Re{\/; i \/, 1 A( 7 )A(z>)*} 



(D.29) 



(D.30) 



x[g£ x - - (q 2 + E y +E x -) cos9 + g£ x + cos 2 # - ffl y +m L cos0], (D.31) 

cos 2 0) 



^-|A(Z)|X 2 [(L 2 c L + i? 2 c,)(|0:f | 2 - \0'i\ 2 )q{E x - + E x 
+ (L 2 c L - J R 2 c i? )2 J Re{o; j L 0;f }m x +m x - costf 
+ (L 2 c L - R 2 e c R )(\0' t f\ 2 + \0\f \ 2 )(q 2 + E xf E x -)cos9 



9 



-L e c L E 2 Re{v* l V jl A(Z)A(i>y[0' l fm x+ m- cos 9 
cos" t>w 1 4 j 

-0; L (gE x7 - (g 2 + E x +E x -)cos^ + gE x + cos 2 #)]}, 

-jc L |^| 2 |^| 2 |A(,>)| 2 E 2 

x [qE - — (q 2 + E +E -) cos + + cos 2 9} . 

A A Aa A A A. A 



(D.32) 



(D.33) 



(D.34) 
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D.2 Chargino decay into sneutrinos 

For chargino two-body decay into sneutrinos 

xt(P x f, Ai) - t + v t \ t = e, fi, r, (D.35) 
the chargino decay matrix is given by 

Pn(xt)K^ = ^ D + °K^D- P.36) 
For the chargino decay into an electron or muon sneutrino the coefficients are 

D = |W| 2 K+-<), (D-37) 

^d= M 9 2 \V tl \ 2 m xt (s a xr p e ); for £ = e, fi, (D.38) 

where the sign in parenthesis holds for the conjugated process x~ £~&£. For the 
decay into the tau sneutrino the coefficients are given by 

D = ^(|ya| 2 + n 2 |^ 2 | 2 )K+ -<), P-39) 

= { +)9 2 (\Va\ 2 - Y T 2 \U i2 \ 2 )m xt (s a xr p T ), (D.40) 

where Y T —m T / (\^2m w cos (3) is the r Yukawa coupling, and the sign in parenthe- 
sis holds for the conjugated process %~ ~^ t~v t . 



D.3 Chargino decay into the W boson 

For the chargino two-body decay into the W boson 

Xt(p x +, \) - Xn(Px°, *n) + W + ( Pw , A fc ), (D.41) 

the decay matrix is given by 

Pih (xt)yx} = E T%%T$f, (D.42) 

* A * 
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with helicity amplitude 

T X D :% = igu(p x o, \ n )r[0 L ni PL + OKp R ]u(p xp \i)e* k *. (D.43) 
For the subsequent decay of the W boson 

W + ( Pw , A fc ) - f'( Pf ,\f>) + f(Pf, A/), (D.44) 
the decay matrix is 

PD 2 (W + ) KXk = £ T D ! J k T D [J k (D.45) 

V ' x f 

and the helicity amplitude 

T D (Z = i^fbf » *f)-fPLv(p f , A/)< fe . (D.46) 

The polarization vectors e^ fc , A fc = 0, ±1, are defined in (B.38). With the set of 
chargino spin vectors s a +, given in (B.14), we obtain for the chargino decay matrix 

pdM)$ = (tyw^r + <^ a D r)< k ^' k (d.47) 

and for the W boson decay matrix 

PD 2 (W + ) KXk = D?e^\ (D.48) 

The expansion coefficients are 

D? = g\\0«\ 2 + \0^\ 2 )[2p^ xt - (jfyrfy + tf x +f w ) - \{m xt + mjo - m 2 w )g^\ 
+2g 2 Re{0**0 L ni )m xt m xl g^ ^ig 2 (\0«f - \0 L J 2 )e^ v x +, a Pw* P-49) 
S" D r = { ±)?(W - \0 L J)m xt [s a x i (p xf -f w ) + s^ip^ -p» w ) + (^ r p w )^] 
-^ 2 (|QS| 2 + \0 L J 2 )m xt e^s xh > x + p - p w , p ) 



+2ig'Re{0«:0^)m x ^s a xta p xt ^ 

a, ix v a.v 
> + p + — 5 

Xj Aj X 



2^ 2 7m(0^0^)m x o«fp^ + - ( eQ123 = 1), (D.50) 
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and 



D? = g 2 (-2pffi + fwf] + Pftw ~ \m 2 w gn M i9 2 e^Pw,« p fJh (D.51) 

where here, and in the following, the signs in parenthesis hold for the charge con- 
jugated processes, xT ~> W~x® an d W~ — > /'/, respectively. In (D.47) and (D.48) 
we use the expansion (E12) for the W polarization vectors 

e^ey = \5^I, V - ^^P P wt C w a (J c ) KXk - ¥wJwA JCd ) KXk - (D.52) 

The decay matrices can be expanded in terms of the spin matrices J c and J cd , 
given in Appendix R2. The first term of the decay matrix p Dl (D.47), which is 
independent of the chargino polarization, then is 

D^e^ k *et' k = IV AfcA ' fc + c £>i(J c ) AfcVfc + cd D 1 (J cd ) XkX 'k , (D.53) 
with 

(m 2 + - mlo) 2 



Di = \g\\0^ \ 2 + \0 L J 2 ) m' x+ + m^o - 2m 2 w + 



m 



w 



cd 



-2g 2 Re(0^0^)m xt m x o, (D.54) 
s £>! = ( y(|0* | 2 - \0 L J 2 )m w {t c w ■ p 4 ), (D.55) 

Di = -g\\o*\ 2 + \o L j 2 ) [(^■p X(+ )(4-p X(+ ) + iK + + 4-^ 

+g 2 Re{0^0 L ni )m xt m xl 5 cd . (D.56) 

As a consequence of the completeness relation (F.14), the diagonal coefficients are 
linearly dependent 

11 D 1 + 22 D 1 + 33 D 1 = -§£>!. (D.57) 

For large chargino momentum p +, the W boson will mainly be emitted into the 
forward direction with respect to p +, i.e. p + pz pw, with p = p/|p|. Therefore, 
for high energies we have (t]^ ■ p +) « in (D.56), and in 11 D 1 ^Z^. 
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For the second term of po 1 (D.47), which depends on the polarization of the decay- 
ing chargino, we obtain 



spa fiu A fc *^ A fe 

with 



spa 
lj D 1 



cspa 



W(|o*l 2 -l^ 12 



(~)3 



■Pw)[ 



o 

m + 



m 2 o 



2m 2 w 



(D.58) 



(D.59) 



(|OSr + l^| 2 K+-2 J Re(0^ 



X 



2 9 2 



cdsr^a 



{ ^)^ 2 (|0^| 2 -|0^| 2 )m x+ x 



b + P Y + P\V L W ' 



(s a xf ■ PwW d - • p xf )(t d w ■ s xf ) - (4 • p xf )(t c w ■ s xf ) 



(D.60) 



(D.61) 



Inserting (D.53) and (D.58) into (D.47), we obtain the expansion of the chargino 
decay matrix in the scalar (first term), vector (second term) and tensor part (third 
term): 



(5 A ^ cd J D 1 + ( 7 A , Ai c ^ i )(J cd )^. 

A similar expansion for the W decay matrix (D.48), results in 

Pd 2 (W + ) x' k x k = D 2 5 X '^ + C D 2 (J c f^ + cd D 2 (r d f^, 

with 

D 2 = \g 2 m 2 w , 
C D 2 = i+^mwitw-Pf), 



cd D 2 



(t c wPf)(t d w -p f )-lm 2 w 5 cd 



(D.62) 

(D.63) 

(D.64) 
(D.65) 

(D.66) 



The diagonal coefficients are linearly dependent 



U D 2 + 22 D 2 + 33 D 2 = -\D 2 . 



(D.67) 



Appendix E 



Neutralino and Chargino two-body 

decay widths 



For the two-body decay of a massive particle in its rest frame 

a-^b + c (E.l) 

the decay width of particle a is 



T(a^bc) = —— — - \T\ dQ — - — — T , (E.2) 

32 ix l m l a J 16 7r m 6 a 

with X(x, y, z) = x 2 + y 2 + z 2 — 2(xy + xz + yz) and \T\ 2 the amplitude squared for 
decay (E.l), where we average over the spins of particle a and sum over the spins 
of particles b,c. 



E.l Neutralino decay widths 

We give the tree-level formulae for the neutralino two-body decay widths T x o for 
the decays 

X°i -»• e RjL e, finLH, f m T, v t v t , x° n Z, xZ w± , xl H h I = e,fi,r; m = 1,2. (E.3) 



E.l Neutralino decay widths 
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Neutralino decay into right selectrons or smuons: xl — > £r + £ ; £ = e,fi 
|T| 2 (X° - l\n = y l/^| 2 Ko - ml), (E.4) 

( TTh^ 777^ 1^ 

r(x° - 40 = f— - ^ ,Uf| 2 . (E.5) 

X 



Neutralino decay into left selectrons or smuons: Xi — > + £ ; £ = e, fi 

\T\\t - ftr) = ^|^| 2 K„ - mfj, (E.6) 

f 777^ 777^ 1^ 



Neutralino decay into staus: X°^t+ + t ; m = 1,2 

m 2 (x° - f+r") = ^(|aLl 2 + lCJ 2 )Ko - m 2 J, (E.8) 
rfe° - f+r") = (m g""f J V (|c4il 2 + l^l 2 )- (E-9) 

OZ 7T III n 

Xi 



Neutralino decay into sneutrinos: Xi —> v n + £ = e,fi,r 

\T\\t - ^) = y |/^| 2 Ko - ml), (E.10) 

(777^ 777^ 1^ 

r(x?^) = |— 5 "S W- (E.n) 

OZ /i //6 
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Neutralino and Chargino two-body decay widths 



Neutralino decay into Z boson: x% ~* Z + x n 



~0> 



cos 2 6w 



6m x om x o[(ReO:f) 2 -(ImO' n tr\ + 



"L\2 



2 , 2 9 2 , K»- m X») 2 

m v o + m v o - 2m z H l - — = 



J\(m 2 o,m 2 , rn 2 7 ) 



16 7r m 3 

x" 



Neutralino decay into W boson: x° — ► + Xj 



+ ^(|OSl 2 + |Ojf) 



(m 2 o - m^ T 



2 ^2 



w 



16 7r m 3 

x- 



T| 2 (x° - W+X- ). 



An 



(E.12) 
(E.13) 



(E.14) 



(E.15) 



Neutralino decay into Higgs boson: xf - ► 

\T\\x° - Hlxl) = 2g 2 m x om x o[Re(H^Re(H^ + Im(H^)Im(H^\ + 
9 2 



\(m 2 x o + m 2 x0n -m 2 H o)(\H^ + \H^), 
./A(mVmLraL) 

1U Ji III n 



(E.16) 
(E.17) 



with H% = Q"* cos a - 5j sin a, Hjj = Hjf and 
1 

13 ~ 2cos9 w 

S" = 1 
13 2cos6 w 



{(N l3 cos + N iA sin (3)N j2 + (iV i3 cos (3 + iV j4 sin /3)iV i2 ] , (E.18) 
[(NH cos /3 - N i3 sin /?) iV i2 + (N j4 cos /3 - N j3 sin /3) iV i2 ] . (E.19) 



E.2 Chargino decay widths 
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The Higgs mixing angle a for small tan (3 can be obtained approximately by 
diagonalization of the neutral Higgs mass matrix 



M 



m 2 z cos 2 (5 + m\ sin 2 f3 
— (m| + m 2 A ) cos f3 sin (3 



[m 



| + m 2 A ) cos /9 sin j3 \ 



m A cos 2 /3 + m| sin 2 (3 + 5 t / 



(E.20) 



which includes the largest term (top-loop) of the one-loop radiative correc- 
tions 



St = 



3g 2 mf 

Ifi 2 2 • 2 / o 1 °g 

iO 7T Tflyy Sin /? 



'm^m 2 ^ 



(E.21) 



We obtain for the Higgs masses 



(m H o) 2 = - |m* + M 2 ^ - v/(M£ - Mg)* + 4 (Mg)* 
1 



+ M* + */(Mg - Ml) 2 + 4 (M^) 2 



(E.22) 
(E.23) 



For the mixing angle we obtain 



cos a = 



sin a 



(MR) 2 + [MR - {m H o) 2 \ 2 

MR ~ (™hQ 2 
^(MR) 2 + [MR - (m„o) 2 ] 2 



(E.24) 
(E.25) 



If we choose a large Higgs mass parameter, e.g. m A = 1 TeV, we have very 
approximately m H o fx 115 — 130 GeV and m H o fx M a , which follows from 
(E.22) and (E.23). In addition, explicit CP violation is not relevant for the 
lightest Higgs state [67]. 



E.2 Chargino decay widths 

We give the tree-level formulae for the chargino two-body decay widths T + for 
the decays 

xt -> W + x° n , el"e, fit^, t£ 2 v t , e + u e , /i + z> M , t + u t . (E.26) 
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Neutralino and Chargino two-body decay widths 



For the heavy chargino xt a ^ so the decays into the lightest neutral Higgs boson H® 
and the Z boson are possible 

Xt - xtZ, xtK (E.27) 
• Chargino decay into W boson: xt —> W + + Xn 

\T\\xt - W+xl) = 9 j(\0«f + \0 L J>) [ml f + mjo - 2m* w + H + " ^ 



w 

-6g 2 m x+ m x oRe(0^0^), (E.28) 



r(xt - W+xZ) = V |T| 2 (x+ - (E.29) 



Chargino decay into selectrons or smuons: xt —> £\ + ^ £ = e,/i 
\T\ 2 (xt - %v e ) = y l^i| 2 K i+ - mf), (E.30) 
(m 2 + - mf) 2 

r(xt - ft^) = ^9 2 \Ua\ 2 - (E.31) 



Chargino decay into staus: xt ~^ T m + v r\ m = 1,2 

\T\\xt - r+i/ T ) = ^|CI 2 K + - <), (E-32) 



r(x+ - f> T ) = g 3 ^ g'ici 2 , (E-33) 

and defined in (A.37). 

Chargino decay into electron or muon sneutrinos: xt —> £ + + &e; £ — e,/j, 

\T\\xt - Put) = yl^i| 2 (< + - m£), (E-34) 
(m 2 + - ml) 2 

^-^--g^-Wrf. (E.35) 



E.2 Chargino decay widths 
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Chargino decay into tau sneutrino: xt —> T+ + v r 

\T\\xt^T+v T ) = 9 4(\Va\ 2 + Y T W)(ml + -ml), 



T V T 



(m 2 + — ml) 2 



and Y T defined in (A.40). 

Chargino decay into Z boson: xt —> Z + Xi 

\T\\xt - Zxt) = -|^-{ - 6 m xtmxt Re(0«M 2 ) 



+ ^(|0f 2 | 2 + \0 L 12 \ 2 ) [m 2 xt + m\ x - 2m| + ^ 



J \(m 2 + , m 2 i , mi) 

nxt - zxt) = v * 2 .J 1 m 2 (^ - z^). 



16 7r m 3 + 

*2 



(E.36) 
(E.37) 



(E.38) 
(E.39) 



Chargino decay into Higgs boson: xt —> + xt 

\T\\xt - Hlxt) = 2g 2 m xt m xt [Re(F 1 L 2 )Re(F 1 R 2 ) + Im{F^)Im{F^)] + 



^K.+< + -< ? )(i^i 2 + ra 2 ), 



16 7r m d 



X2 



with itf = ^(E/£V£ sin a - [/* J$ cos a) and = 



(E.40) 



A /A(m 2 + ,m 2 + ,m 2 ) 
nxt - fl?xJ") = V X 3 1 l T ! 2 (^ - Hlxt\ (E.41) 



The Higgs mixing angle a is given in (E.24) and (E.25). 



Appendix F 



Spin formalism for f ermions and 

bosons 



F.l Bouchiat-Michel formulae for spin | particles 

For the calculation of cross sections we expand the spin density matrices in terms 
of the Pauli matrices, see e.g. (C.l), (C.10) for neutralinos. This expansion is straight 
forward if for the neutralinos or charginos a set of spin-basis vectors s a ^ has been 
introduced, see (B.14). Together with j3 M = p^/m they form an orthonormal set 



p ■ s a = 0, (F.l) 

s a - s b = -5 ab , (F.2) 

*X = -9r + ( sum over «)■ (F-3) 
The helicity spinors are normalized by 

u(p,\)u(p,\') = 2m5 xy , (F.4) 

v(p,X)v(p,X') = -2m5 xx , (F.5) 

The Bouchiat-Michel formulae for massive spin 1/2 particles are then [68] 

u(p, A') u(p, A) = \ [5\\r + 75 fo^W + m), (F.6) 

v(p, A') v(p, A) = ^ [5 A 'a + 75 ^a'aK^ - ™)> ( sum over «)■ ( R7 ) 



F.2 Spin formulae for spin 1 particles 
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F.2 Spin formulae for spin 1 particles 



The Bouchiat-Michel formulae for spin 1/2 particles can be generalized for higher 
spins [69]. In order to describe the polarization states of a spin 1 boson, we have 
introduced a set of spin vectors t a , see (B.37). Note that they are are not helicity 
eigenstates like the polarization vectors e* k , defined in (B.38). The spin vectors 

and k^ = k^/m form an orthonormal set: 



k ■ t a = 0, (F.8) 
t a . t b = -5 ab , (F.9) 

tl K = + JL i; (sum over a). (RIO) 

7TZ 



The 3x3 spin 1 matrices J c obey [ J c , J d ] = ie c d e J e and are given below. We can 
define six further matrices 



jcd = jcjd + jd jc _ ^cd (R11) 

with J 11 + J 22 + J 33 = 0. They are the components of a symmetric, traceless tensor. 
We now can expand [69] 

e*eP = \5 X '^I» V - J- e ^p p z t c >"(J c f^ - \t c / v {J cd ) x ' k x \ (F.12) 
summed over c, d. The tensor 



V = ~9,u + ^ (K13) 



guarantees the completeness relation of the polarization vectors 



The second term of (E12) describes the vector polarization and the third term de- 
scribes the tensor polarization of the boson. 
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Spin formalism for fermions and bosons 



In the linear basis (B.37) the spin-1 matrices are defined as ( J£) jfe = —ie 



cjk- 




T 2 



J? 



t23 



The matrices J c and J cd in the circular basis, see (B.38), are obtained by the unitary 




7 3 



J? 



Tl3 



/0 -i 





(E15) 
(E16) 
(F.17) 



transformations J c 



i 

v/2 



-_L. n % — 

1 



At. J£ 



At 



A and J' 



At • J£ d 



A, respectively with 



i-i 



(E18) 



J 1 



J 



ii 




J 



12 



/ i\ 


V -i J 



J 2 



22 



J 



23 






( 


-i 


0\ 












J 


K 





l) 



J 



33 



J 



13 




(E19) 



(E20) 



(E21) 



In calculating products of density matrices, the following relations are helpful 



Tr{J a } = 0, Tr{J afe } = 0, Tr{ J a J hc } = 0, 



ie 



abc 



Tr{JV 6 } = 25 a \ Tr{ J a J b J c } 
Tr{J ab J cd } = -U ab 5 cd + 25 ad 5 bc + 25 ac 5 l 



(E22) 

Tr{ J a J b J c J d } = 5 ab 5 cd + 5 ad 5 bc , (E23) 

(E24) 



Appendix G 



Definitions and conventions 



We use natural units c = 1, h/2n = 1. 
The metric tensor 



w = (an ■= 



(i o o o \ 

0-10 
0-10 
\0 -I J 



defines scalar products 

(a . b) := g^b v = a"b^ = a°b° - ab 

between covariant and contravariant four-vectors 

a M := (a , a 1 , a 2 , a 3 ) = (a , a), a M := g^a^ = (a , a±, a 2 , 03) = (ao, -a). 

The total antisymmetric e-tensor is defined as 



-fit/per 



+1 , if pupa is an even permutation of 0123, 
— 1, if [iv po is an odd permutation, 
, if any two indices are the same. 



The analog definition in three dimensions, the Levi-Cevita-Tensor is 

+1 , if ij k is an even permutation of 123, 
( ; _ e ijk ._ ^ _i ^ if ij A; is an odd permutation, 

, if any two indices are the same. 



(G.l) 



(G.2) 



(G.3) 



(G.4) 



(G.5) 
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Definitions and conventions 



Useful relations of the e-tensor 

-e a ^e a( s p(7 = 2(5^ - 5*8 v p ), -e a ^e a ^ p = 6 5;, -e a ^e a ^ u = 24, (G.6) 

with <f£ = g^g vp . 
The Pauli matrices are 

CTl = ( ? o ) ' a2= (°7)' ^ 3 = ( o - 1 ) • (G - 7) 

The Dirac matrices obey the commutation relations 

{Y, Y} = YY + 7 V = 20"". (G.8) 
In the Dirac representation they read 

and 75 := -?7o7i7273 = ?7°7 1 7 2 7 3 = 7 5 . 
Trace theorems: 

PL,fl:=|(l=F75), ^:=7%, Tr{P Lii? } = 2, 

Tr{^P LiiJ } = 0, Tr{^P L)i? } = 2(a-6), Tr{^ y ^ P L , R } = 0, (G.10) 

Tr{^ y <?$Pl,r} = 2 [(a-b)(c-d) - (a-c)(6-d) + (a • d)(b ■ c)] =p 2i[a, b, c, d], 

and [a, b, c,d] := e fiupa a fl b'' 'c p d° '. 

For numerical calculations we have used the values 

a = 1/128 fine — structure constant at 500 GeV 

sin 2 6 W = 0.2315 weak mixing angle 
mw = 80.41 GeV W boson mass (r -\-\\ 

T w =2.12 GeV W boson width { ' 

mz = 91.187 GeV Z boson mass 
Y z = 2.49 GeV Z boson width 
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